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Abstract 

Aims. The HIFI instrument onboard Herschel has allowed high spectral resolution and sensitive observations of ground-state transi- 
tions of three molecular ions: the methylidyne cation CH^, its isotopologue ''CH^, and sulfanylium SH*. Because of their unique 
chemical properties, a comparative analysis of these cations provides essential clues to the link between the chemistry and dynamics 
of the diffuse interstellar medium. 

Methods. The CH^, '^^CH^, and SH* lines are observed in absorption towards the distant high-mass star-forming regions (SFRs) 
DR21(0H), G34.3+0.1, W31C, W33A, W49N, and W51, and towards two sources close to the Galactic centre, SgrB2(N) and 
SgrA*+50. All sight lines sample the diffuse interstellar matter along pathlengths of several kiloparsecs across the Galactic Plane. 
In order to compare the velocity structure of each species, the observed line profiles were deconvolved from the hyperfine structure 
of the SH^ transition and the CH+, '^^CH^, and SH^ spectra were independently decomposed into Gaussian velocity components. To 
analyse the chemical composition of the foreground gas, all spectra were divided, in a second step, into velocity intervals over which 
the CH+, "CH^, and SII+ column densities and abundances were derived. 

Results. SH^ is detected along all observed lines of sight, with a velocity structure close to that of CH* and '^CH*. The linewidth 
distributions of the CH+, SH+, and '^CH+ Gaussian components are found to be similar. These distributions have the same mean 
((Af) ~ 4.2 km s"') and standard deviation {<r{Av) ~ 1.5 km s"'). This mean value is also close to that of the linewidth distribution 
of the CH^ visible transitions detected in the solar neighbourhood. We show that the lack of absorption components narrower than 
2 km s"' is not an artefact caused by noise: the CH^, '^^CH^, and SH^ line profiles are therefore statistically broader than those of 
most species detected in absorption in diffuse interstellar gas (e. g. HCO^, CH, or CN). The SH+/CH+ column density ratio observed 
in the components located away from the Galactic centre spans two orders of magnitude and correlates with the CH^ abundance. 
Conversely, the ratio observed in the components close to the Galactic centre varies over less than one order of magnitude with no 
apparent correlation with the CH^ abundance. The observed dynamical and chemical properties of SH^ and CH^ are proposed to trace 
the ubiquitous process of turbulent dissipation, in shocks or shears, in the diffuse ISM and the specific environment of the Galactic 
centre regions. 
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1. Introduction 

Studying the diifuse phases of the interstellar medium (ISM) 
is essential, not only because they contain a large part of 
the total gas mass of the cold ISM and are the precursors of 
dense clouds, but also because they harbour the first steps of 
interstellar chemistry. Since the detection of the first diatomic 
molecules CN, C!H, an d CH^ (see references in the review of 
ISnow & McCal]||2006h through their narrow visible absorption 
lines, the improvement of the observational techniques and 
instruments has provided a more comprehensive view of the 



* Based on observations obtained with the HIFI instrument onboard 
the Herschel space telescope in the framework of the key programmes 
PRISMAS and HEXOS. 

** Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



diffuse ISM and led to a better understanding of its dynamical, 
thermal, and chemical evolution. A wide variety of diatomic 
and triatomic molecular species has akeady been observed 
in the diffuse medium. Moreover, its chemical composition 
has now be en probed in the solar n eighbourhood through 
UV (e.g. IShull&BeckwithI 119821: ICrawford&WiUiarnsI 



19971 ISnow et al J 120001: iRachford et al.1 l2002t iGrv et alJ l2002t 



1995'; 'Gredell 
Maier et al 



Lacour et al. visible (e.g. Cran e et al 

1997; Thorburn et al. 2003; Weselak et all 
2001), and radio (e.g. |Haud& Kalberli 
i2006 and references therein) spectroscopy, and the inner 
Galaxy material through submillimet re and radi o (millime 
tre, centimetre) spectroscopy (e.g. Koo 1997; Fish et al 



20081: ^ 

120071; [Liszt et al 



i 



2003^__>Ivman & Mijlar 1989; Greav es & Williams. ^99; 
iNeufeld et al.ll2002HPlume et alJl2004 . 



The Herschel space mission has broadened this investiga- 
tion, giving access to the full submillimetre domain, which has 
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allowed the detection of many molecular species that could not 
be detected from the ground before because of the high opacity 
of the atmosphere. In the framework of the HIFI key programme 
PRISMAS (PRobing InterStellar Molec ules with Absorption 
line Studies) many hydrides such as HF (Neufeld et al.ll 20l"0bl 
Sonnentrucker et al. 2010), OH+, Ha O^dNeufeld et al. 201Qi 
Gerinetal. 2010a) CH (jGerinetal. 2010b), NH, NHj, NH3 
(Persson et al. 201(1), and CH+ (Falgarone et al. 2010) were de- 
tected in absorption against the strong continuum emission of 
distant star-forming regions, providing for the first time a good 
census of these light hydrides in the inner Galaxy. 

Of all molecules targeted by PRISMAS, the methylidyne 
cation CH^ is particularly interesting because its presence in 
the diffuse ISM remains one of the most intriguing puzzle in 
astrophysics. The CH+ abundances predicted by steady-state, 
UV-dominated, PDR-type (PhotoDissociation Regions) chemi- 
cal models are several orders of magnitude lower than the ob- 
served values, because all the CH^ formation pathways that are 
alternative to the highly endothermic H- H2 CH^ + H re- 
action, are ineffici ent in balancing its fast destruction by hydro- 
genation. Indriolo ^t al.l (12010.) recently showed that the upper 
limit on the CH^/CH^ abundance ratio observed towards Cyg 
OB2 can only be reproduced in diff'use molecular clouds with 
extreme conditions (i.e. /h, < 0.2, or T > 1000 K ). So far, 
the solution to this puzzle has been to invoke regions of the dif- 
fuse gas where a warm chemistry is activated by the release of 
supra-thermal energy in the cold ISM in duced by low-velocity 
magnetohydrodynamic (MHD) shocks dPrain e & Katz '1986^, 
f'ineau des Forets et al. 1986), Alfven waves (Federman et al. 
[1996). turbulent mixing (Xie et al. 1995; Lesaff're et al. 2007), or 
turbulent dissipatio n (Falgarone et al. 1995: Joulain et al...l998; 
iGodard et aLll2009t) . Undriolo et al. (2010) claimed that observa- 
tions of the excited levels of CH3 are able to provide the clues 
necessary to favour one theory over the other. 

A potentially related problem is the existence of sulfanylium 
(SH^) in the diff'use gas. Because the hydrogenation reaction of 
S^ has an endothermicity twice as high as that of C^, a mea- 
surement of the SH^/CH^ abundance ratio should provide valu- 
able insights about the regions where CH^ is formed. Sought for 
without success in the local diffuse ISM in ab sorption against 
nearby stars in the UV domain since 1988 (Mil lar &Hobbsl 
[1988; Magnani & Salze r 1989, 1991), t his molecular ion was re- 
cently detected by Me nten et alJ (|201 ih through its ground-state 
rotational transition in the submillimetre range observed in ab- 
sorption towards the Galactic centre line of sight SgrB2(M) with 
the APEX telescope. 

In this paper, we report the detection of SH^, CH^, and 
'^CH^ towards the six distant star-forming regions DR21(OH), 
G34.3+0.1, W31C, W33A, W49N, and W51, and the Galactic 
centre sight line^ SgrA*+50 and SgrB2(N), and we perform a 
cross analysis of the observed properties of those three species. 
The observations, obtained in the framework of the key pro- 
grammes PRISMAS and HEXOS, are presented in Sect.E] The 
methods used for the analysis and the results obtained are shown 
in Sects.[3]and|4l respectively. We conclude this work in Sect.|5] 
with a discussion on the chemical and dynamical properties of 
the gas seen in absorption. The comparison with the model pre- 
dictions will be the subject of a forthcoming paper (Godard et 
al., in prep.). 



2. Observations and data reduction 

2. 1 . Observing conditions 
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' The SgrA*-l-50 sight line (GCM-0.02-0.07) corresponds to the 50 
km s"' cloud located in the vicinity of SsrA*, w hich is known to be a 



bright submillimetre source dPowell et 



^SgrA^, w 
alJl999h . 



Figure 1. Original SH"^ spectra (double-sideband antenna tem- 
perature Ta before removing the standing waves) observed to- 
wards SgrA*-H50, W31C, W33A, and W49N in the horizontal 
polarization and for the three different LO settings (in black, red, 
and blue). For more clarity the red and blue curves have been 
shifted from the black curve by 0.2 and 0.4 K (for the SgrA*-H50 
data), and by 0.1 and 0.2 K (for the W31C, W33A, and W49N 
data). 



The observations were carried out from March to October 
2010 and in April 2011 towards the eight submillimetre back- 
ground continuum sources listed in Table [T] (with their Galactic 
coordinates, their distance from the Sun, and their measured 
single-sideband continuum temperature Tc in K at ~ 830 and 
~ 530 GHz). Using the dual-beam switch (DBS) mode (with 
a throw at 3' from the source) and the wide-band spe ctrometer 
(WBS) of Herschel/HIFI (see iRoelfsema et al.ll2012l for a de- 
tailed description of the properties and performances of HIFI), 
we observed 

• the 7 = 1 <- absorption Hnes of CH+ and '^CH+, in the 
upper and lower sidebands of band 3a, 

• and the F = 3/2 ^ 1/2, 5/2 <- 3/2, and 3/2 <- 3/2 hyper- 
fine components of the A^, 7 = 1 , 2 <— 0, 1 absorption line of 
SH^, in the lower sideband of band la. 

The data obtained towards SgrB2(N) were part of the 
full HIFI spectral scan performed by the HEXOS key pro- 
gramme; the corresponding double-sideband spectra were de- 
convolved into single-sideband spectra, including the continuum 
(IComito & SchiUce ,2002.) . Towards the other sources, the ob- 
servations were performed in the framework of the PRISMAS 
key programme; to separate spectral features from upper and 
lower sidebands of the WBS spectrometer, each transition was 
observed using three slightly different settings of the local oscil- 
lator (LO) frequency, adapted to induce a relative velocity shift 
of ~ 30 km s"' between the two sidebands. The spectroscopic 
parameters of the observed lines are listed in Table|2] along with 
the rms noise levels relative to the single-sideband continuum 
intensities obtained with an on-source integration time ranging 
from 1 to 20 min. In these frequency ranges, the WBS resolution 
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Source 


RA(J2000) 


Dec (J2000) 


1 


b 


D" 


Rg 




7'.(K)'' 




(h) (m) (s) 


D (') (") 


n 


n 


(kpc) 


(kpc) 


(830 GHz) 


(530 GHz) 


SgrA*+50 


17 45 50.2 


-28 59 53 


359.97 


-0.07 


8.4 


0.1 


0.4 


0.2 


SgrB2(N) 


17 47 19.9 


-28 22 18 


0.68 


-0.03 


8.4 


0.1 


8.0 


2.9 


G34.3+0.1 


18 53 18.7 


01 14 58 


34.26 


0.15 


3.8 


5.8 


2.3 


0.7 


W31C 


18 10 28.7 


-19 55 50 


10.62 


-0.38 


4.8 


3.9 


1.9 


0.5 


W33A 


18 14 39.4 


-17 52 00 


12.91 


-0.26 


4.0 


4.7 


0.6 


0.2 


W49N 


19 10 13.2 


+09 06 12 


43.17 


+0.01 


11.5 


7.9 


2.4 


0.8 


W51 


19 23 43.9 


+ 14 30 31 


49.49 


-0.39 


7.0 


6.6 


2.7 


1.0 


DR21(0H) 


20 39 01.0 


42 22 48 


81.72 


0.57 


1.0 


8.4 


1.3 


0.4 



" For G 34. 3 +0.1, W31C, W33A , and W49N, the source distances D, and the subsequent Galactocentric distances Rc, are taken from lFish et al.l 

(2003) and lPandian et al.l (12008). who resolved the kinematic distance ambiguity. Distance uncertainties are of the order of 0.5 kpc. 
* Uncertainties on the continnuum temperatures 7^ are of the order of 10% (iRoelfsema et al..,2012i) . 



Table 2. CH+X'E^, '^CH^X'S^, and SH^ X^I, spectroscopic parameters for the observed pure rotational transitions. Numbers in 
parenthesis are power of 10. 



Transition 




g« 


vo 


A 




1/ 


















(GHz) 


(s-i) 


(km s-') 




SgrA*+50 SgrB2(N) DR21(0H) 


G34 W31C W33A W49N W51 


CH+ 


1 - 
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835.137504 


5.83 (-3) 






1.1 (-2) 


9.8 (-3) 


1.5 (-2) 


1.9 (-2) 5.2 (-2) 2.0 (-2) 1.7 (-2) 1.7 (-2) 


"CH+ 


1,1/2- 


0,1/2 


2 


830.215004 


5.83 (-3) 


0.59 


0.50 










'3CH+ 


1,3/2- 


0,1/2 


4 


830.216640 


5.83 (-3) 


0.00 


1.00 


5.5 (-2) 


5.5 (-3) 


1.5 (-2) 


2.0 (-2) 1.5 (-2) 2.9 (-2) 1.4 (-2) 2.0 (-2) 


SH+ 


1,2,3/2- 


0,1,1/2 


4 


526.038722 


7.99 (-4) 


+5.26 


0.56 










SH+ 


1,2,5/2- 


0,1,3/2 


6 


526.047947 


9.59 (-4) 


0.00 


1.00 


5.0 (-2) 


1.5 (-2) 


2.9 (-2) 


2.5 (-2) 9.0 (-3) 2.4 (-2) 1.1 (-2) 1.5 (-2) 


SH+ 


1,2,3/2 - 


0,1,3/2 


4 


526.124976 


1.60 (-4) 


-43.93 


0.11 











" CH+ '^CH+ and SH+ line frequencies are from lAmanol ( 120101) and lSavage et al.l d2004l) . respectively. 

* Velocity shifts of the hyperfine components of the SH^ (1,2 <— 0, 1) line, relative to the F = 5/2 «— 3/2 transition, the reference hyperfine 
transition in the text. 

Intensities of the hyperfine components computed (at Tex — » <») as g„A„;/g„oA„/o, relative to a chosen reference hyperfine transition (desig- 
nated by the index in the previous formula). 

CTi/Tc is the rms noise level divided by the continuum intensities of the spectra at the resolution of 1.1 MHz. 



of 1 . 1 MHz corresponds to velocity resolutions of ~ 0.36 km s 
for the CH^ and '^CH^ transitions, and of ~ 0.57 km s ' for the 
SH+ transitions, and the Herschel HPBW is 26" at 835 GHz and 
41" at 526 GHz. 

The data were cal ibrated with hot and cold blackbodies 
(IRoelfsema et al.l [20121) . reduced using the standard Herschel 
pipeline to Level 2, and subsequently analysed using the 
Herschel Interactive Processing Environment (HIPE v5.1, 
[Ott et al. 2010). The final analysis was pe rformed with the 
G1LDAS-CLASS90 softwar^l (fPetvl l20()5h . and a set of 
Fortran95 numerical routines that we developed. While the sig- 
nals measured in the two orthogonal polarizations that were ob- 
tained with the three LO settings agreed excellently for both the 
CH^ and '^CH^ line observations, the SH^ spectra, displayed 
in Fig.[T] exhibit standing waves (SW) of identified origirQ and 
removed using the HIPE sine wave fitting task FitHifiFringe 



^ See http://herschel.es ac.esa.int/HIPE_dovynload. shtml| for more in- 
formation about HIPE. 

^ See http://www.iram.fr/IRAMFR/GILDASI for more information 
about GILDAS softwares. 

* Most of the obs erved SW ha ve a period of ~ 90 - 100 MHz, iden- 
tified by iRoelfsem a et al.l ( 120 1 2|) as reflections occuring between the 
mixer focus and the cold and hot black bodies. For two spectra ob- 
served towards W49N we had to remove an additional standing wave 
with a period ~ 150 MHz. 



(FHF). Since the inferred detected opacities of SH^ are low, and 
since the period and the amplitude of the waves are similar to 
the size and the depth of the SH^ absorption features (see Fig. 
[TJ, the resulting spectra were sensitive to the FHF input param- 
eters: several plausible solutions were found depending on the 
number of sub-bands taken into account for the fit, the number 
of sine waves to remove, and the use of a frequency mask. We 
estimated a maximal error of 50 % on the deepest absorption 
features of G34.3+0.1 and W51, 30 % on those of W33A, and 
W31C and 10 % on the others. The latter value is comparable 
to that due to the u ncertainties on the bea m efficiency and the 
sideband gain ratio ( IRoelfsema et alj|2012h . 

For each transition and both polarizations, we obtained an 
average spectrum by combining the data from the three obser- 
vations with different LO settings. Because we are interested in 
the velocity structure and the properties of the absorbing gas, the 
spectra in both polarizations were normalized to their respective 
continuum temperature. As in Falgarone et al. (2010), we used 
the saturated shape of the CH+ absorption line profiles to mea- 
sure the sideband gain ratios R at 835.1375 GHz, defined as the 
ratio of the continuum temperatures in the lower and upper side- 
bands. For all spectra with saturated absorption lines, we found 
R ~ 0.9 - 1 and R ~ 0.7 - 0.8 in the horizontal and vertical polar- 
izations, respectively. We finally used these values at 835 GHz, 
and = 1 at 830 GHz and 526 GHz, to combine the data from 
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both polarizations, and obtain the final average spectra (with rms 
noise levels given in Table |2]l shown in Fig. |2] as functions of 
the LSR (local standard of rest) velocity. This figure illustrates 
the quality of the baselines over the large bandwidth for most 
of the spectra. The strong emission lines detected in the spectra 
shown in Figs.[r|and|2]were identified as the HCO^ (6 — » 5) and 
H2CO (6 5) transitions near 535. 1 GHz and 525.6 GHz, three 
methanol Hnes near 829.9 GHz and 830.3 GHz, and the SO2 
emission band near 835 GHz, emitted by the SFRs themselves. 



2.2. Deconvolution of the SH+ hyperfine structure 

Unfortunately (see Table |2]i, the velocity shifts Au/, associated 
with the SH^ hyperfine transitions are smaller than the observed 
velocity ranges of the SH^ absorption spectra, and prevent us 
from performing a direct cross comparison of the velocity pro- 
files of the SH+, CH+, and '^CH+ lines. To solve this problem, 
we developed a numerical procedure to extract the signal associ- 
ated with each hyperfine transition, solving the following set of 
equations for Triv), over the entire absorption velocity domain. 



2 Ih(k) Tr(v - Avh(k)) = -In [T(v)/T,] ., 



(1) 



k=l 



where Nh, //,, r,., and T(u)/Ti- are the number of hyperfine tran- 
sitions, the intensities relative to the strongest hyperfine compo- 
nent, the opacity of the reference hyperfine transition, and the 
normalized line profile (line/continuum), respectively. The re- 
sulting spectra are shown in Fig. [3] and, as an example, the out- 
come of the hyperfine decomposition code, applied to the ab- 
sorption lines observed towards SgrB2(N) and W49N, is shown 
in Fig.m This figure illustrates the excellent agreement between 
the original data (in black) and the spectra rebuilt after decom- 
position (in green). 

TheJfCH_^_7 = 1 <— line also exhibits a spin-rotation split- 
ting dAmandbOlOl) . although the associated F - 3/2 «— 1/2 
and F - 1/2 <— 1/2 transitions are separated by only 1.636 
MHz (~ 0.59 km s ') and are therefore too close to be individu- 
ally resolved given the significant velocity dispersion of the gas. 
This hyperfine structure induces a systematic broadening of the 
absorption velocity components that depends on their FWHM 
Aureal : for Afreai Varying between 2 and 10 km s"' the broaden- 
ing ranged between 4% and 0.1%. Since this error is far smaller 
than that imputable to the rms noise levels of the '^CH^ spec- 
tra (see Sect. [3]), we chose to ignore the hyperfine structure of 
'^CH^ in the following analysis. 

The spectra shown in Fig. [3] are highly structured and have 
the following remarkable properties: (1) thanks to the high sen- 
sitivity of the HIFI receiver, the SH^ ion is seen in absorption 
along every line of sight; (2) all hydride lines are detected in ab- 
sorption, and within the limits imposed by the signal/noise (S/N) 
ratio, CH^ '^CH^ and SH^ absorptions are detected over the 
whole velocity range of the foreground matter along each line of 
sight; (3) although the opacity ratios vary from one line of sight 
to another and from one velocity range to another, the velocity 
structure of SH^ is similar to those of CH^ and '"'CH^. It is the 
similarity and differences of these absorption line profiles that 
are the focus of the present study. 



' This result on the line profile broadening is derived from the analy- 
sis of 1760 synthetic spectra taking into accou nt the hyperfin e structure 



3. Analysis of the line profiles 

3.1. Multi-Gaussian decomposition 

As in iGodard et al.l (I2OIOI) . the decomposition of the spectra 
in velocity components was deduced through a multi-Gaussian 
fitting procedure that we developed, based on the Levenberg- 
Marquardt algorithm, which takes advantage of the information 
carried by the hyperfine structure of a given transition. This algo- 
rithm is applied to adjust, in the least-squares sense, the minimal 
number of Gaussians required to describe the data within the ob- 
servational errors without introducing any systematic effect. The 
number of Gaussians was increased, for instance, when we vi- 
sually spotted serpentine curves - characteristic of poor fits in 
the line wings - in the residuals. Thus, for each transition, the 
observed normalized line profile (line/continuum) is written 



T{v) 



- exp 



N,, 

■Y,Yj^h(k)Ti,{j)e 

j=i k=\ 



i.-i.Ql./)-A..,,(H 

I 'o^n I 



(2) 



of '^CH^ (line strength and velocity structure. I Amam 



ti ypernn e 



where tq, vq, and ctq are the usual Gaussian fit parameters. All 
spectra were decomposed independently from one another, with- 
out imposing any constraints on the Gaussian parameters. The 
choice of the input parameters, namely the initial values of voi j) 
and o-qH) for each velocity components j, was guided by the 
comparison of the different lines observed towards each source. 
To correctly determine the opacity of weak absorption features 
blended with saturated lines, as observed in the CH^ spectra, we 
applied an empirical model constrained by the wings of the satu- 
rated line profiles. The results of the multi-Gaussian decomposi- 
tion procedure and the associated errors on the Gaussian param- 
eters are given in Table lA.ll of AppendixlAl and the resulting fits 
and models of the saturated line profiles are displayed in Figs. 
IA.1I - IA.8] Because we aim to compare the kinematic signatures 
of the CH^, '^CH^, and SH^ spectra, we discuss below the relia- 
bility of the extracted Gaussian components in view of the errors 
on their parameters. 



3.2. Validity and self-consistency of the multi-Gaussian 
decompositions 

Since the decomposition algorithm allows the detection of com- 
ponents with very weak central opacities, the numerical proce- 
dure may converge upon Gaussian components whose reality is 
questionable. To keep only the most reliable velocity compo- 
nents for our subsequent analysis of the linewidths, we applied 
the following detection criterion: any Gaussian component was 
considered real if its Gaussian parameters simultaneously ver- 
ify Au > 3cr{Av) and t > 2.5cr(T), where o-(Av) and o-(t) are 
the errors on Av and r respectively. The resulting confirmed or 
uncertain Gaussian components (in the following C-components 
and U-components) are indicated in Table lA.ll dn columns 4, 8, 
and 12), and in Figs. lA.l l- IA.8l (solid red and dashed blue curves) 
of AppendixlAl 

In total we found that 25 C-components are simultaneously 
observed in at least two molecular spectra towards DR21(OH), 
G34.3+0.1, W31C, W33A, W49N, and W51. When compared, 
the positions of these components are found to agree with one 
another within their respective error for 18 of them, within 0.5 
km s ' for 6 of them, and within 1.3 km s ' for 1 of them; sim- 
ilarly, out of the 17 common C-components observed towards 
SgrA*+50 and SgrB2(N), 10 are found to agree with one an- 
other within 1 km s"', 6 within 2 km s"', and 1 within 3.5 km 
s"' . Except for the SH^ components at -126 km s"' observed to- 
wards SgrA*-i-50 these shifts in the central positions are at least 
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four times smaller than the corresponding Gaussian linewidths. 
Moreover, many of the U-components have corresponding C- 
components in other species, indicating that the fitting process is 
robust and that the selection method is severe. 

These concordances combined with the strict selection on the 
Gaussian parameters suggest that all C-components are real de- 
tections and not artefacts caused by noise or the standing waves 
removing procedure. 

3.3. Comparison of the Gaussian iinewidtlis 

In Fig. |5] we display the distributions of linewidths associated 
with the C-components extracted from the CH^, '^CH^, and 
SH"^ spectra. To emphasize the differences observed along the 
Galactic centre sight lines (/ ~ 0), these distributions are com- 
puted for all sight lines (left panels), and for the / sight 
lines only (right panels). As a comparison, the linewidth distri- 
butions of the HCO ^ ground state radio transition observed by 
iGodard et al.1 jloTol) towards W31C, W49N, W51, and G34.6 
are displayed in panels (a), and those of CH and CH^ visible 
transitions observed at high sp ectral resolution (~ 0.3 km s"') 
in the local diffuse medium bv lCrane et alj ( Il995h are shown in 
panels (c) of Fig.|5] Lastly, panels (d) and (e) of Fig.|5]display the 
first- and second-order moments of the Au-distributions issued 
from the combined CH^, '^CH^, and SH^ data (red squares) and 
from the CH^ data observed in the solar neighbourhood (blue 
circles). In panels (a), the histogram of HCO^ linewidths is nar- 
rower and peaks at lower values than those of the submillime- 
tre lines of CH^, '^CH^, and SH^. Conversely, the histogram of 
the visible CH^ data, characterising the local diffuse interstel- 
lar matter, is very similar to that of the submillimetre data. We 
discuss the validity and the significance of these comparisons 
below. 

To demonstrate that the absence of narrow velocity compo- 
nents in the CH^, '^CH^, and SH^spectra is real and not due to a 
limitation of our extraction algorithm, we have derived the mini- 
mum width of a Gaussian component of optical depth tq that can 
be extracted from a given profile characterised by a noise ctj at 
the velocity resolution 6v. Because ctj scales inversely with the 
square root of the velocity resolution, a component detected at a 
3cr^ level necessarily verifies 

Av > 5v ['iarlTof . (3) 

While the noise level of the SH^ spectra towards SgrA*-H50, 
DR21(OH), G34.3+0.1, and W33A forbids us to extract compo- 
nents with linewidth smaller than 1 km s ', that of all the other 
spectra is small enough to allow us to detect narrow velocity 
structures down to two velocity channels (0.72 km s"' for CH^ 
and '''CH^, and 1.14 km s ' for SH^). Therefore the scarcity of 
components in the first bin [0 ; 2 km s '] of the CH^, '^CH^, 
and SH^ histograms is not a noise artefact. 

Now, can we compare these linewidth distributions with 
one another? Beca use the S/N ratios of HCO^ (see Table. 1 of 
iGodard et al.ll2010l) and CH^ (see Table. |2]i profiles are high and 
comparable, those two sets of spectra are decomposed with the 
same level of detection. The noise therefore affects the statistics 
of the component linewidths at the same level. The same is true 
for the lines of SH^ and '^CH^, which have comparable, though 
poorer, S/N ratios. Finally, because CH^ and '^^CH^ necessarily 
bear the same dynamical signatures, all distributions displayed 
in Fig.|5]can be compared with one another 

We recall that all absorption spectra have been decomposed 
into Gaussians independently from one another, without impos- 
ing the same velocity centroids or width to the different velocity 



components of a given line of sight. In addition, since the most 
intense SH^ components are saturated in CH^, and because the 
moderate S/N ratio of the '^CH^ and SH^ data prevents us from 
observing the low-opacity structures detected in the CH^ spec- 
tra, the respective distributions do not correspond to the same 
velocity components. It is therefore remarkable that the distribu- 
tions of the CH^, '^CH^, and SH^ widths of the Gaussian veloc- 
ity components are so similar, even identical within the statistical 
uncertainty. They exhibit a common pattern: a peak (hereafter 
called PI), observed towards all the sources, defined by a first- 
(mean) and a second- (standard deviation) order moment^ of 
4.2 + 0.2 km s"' and 1.5 ±0.1 km s"', respectively, and an ex- 
tended tail (hereafter called P2), observed only on the Galactic 
centre sight lines (left panels), with Aw-values up to 20 km s"' . 

Because the lines of sight sample kiloparsecs of interstellar 
material in the Galaxy, the patterns PI and P2 result from the 
small-scale dynamics of the production processes of CH^ and 
SH^, the turbulent dynamics of the diffuse ISM, and the Galactic 
dynamics. The comparison between the left and right (d) and (e) 
panels of Fig. |5] shows that the first- and second-order moments 
of the Aw distributions are the same for the components of the 
Galactic ISM along the Z sight lines and the visible CH+ 
lines sampling the solar neighbourhood (defined by first- and 
second-order moments of 4.3 ± 0.4 km s ' and 1.85 + 0.3 km 
s"' , respectively). Because the latter is unaffected by the Galactic 
dynamics, this similarity suggests that the peak PI results from 
the dynamics of the formation processes of CH^ and SH^ con- 
volved with that of the turbulence of the diffuse gas. Because the 
tail of the Ai/ distributions is observed only on the Galactic cen- 
tre sight lines, and because the high absorption dips (Aw > 8 km 
s"') observed on the / sight lines can b e decomposed into 
many narrow components (as performed by IGodard et alj|2010l 
with the HCO^ spectra), we conclude that P2 is caused by the 
Galactic dynamics only. 

Lastly we note that while the first- and second-order mo- 
ments of pe ak PI both agree w ith those of the Aw distribution 
obtained by ICrane et aP (Il995h within the statistical uncertain- 
ties, they also substantially differ from those o f the CH^ Af dis - 
tribution derived in the solar neighbourhood bv lPan et"an (l2004l) . 
who find first- and second-order mome nts of 3.3 + 0.04 km s ' 
and 0.4 ± 0.03 km s As proposed bv lPan et all (12005,) . these 
differences may originate from their profile fitting method: (1) 
the number of CH^ components and their position were con- 
strained by the observations of other species, possibly unrelated 
chemically (e.g. CN, CO, Ca), and (2) a maximal value of Ay of 
5.8 km s"' was set. 

4. Analysis of the integrated opacities 

Since the multi-Gaussian decomposition of CH^, '^CH^, and 
SH^ absorption profiles provides sets of C-components that do 
not always strictly coincide in velocity and width, and sets of 
U-components that sometimes clearly correspond to real absorp- 
tion but with width and depth poorly constrained, our subsequent 
analysis and comparison of the column densities and abundances 
of these species is based on integrals of opacities computed over 
broad velocity intervals corresponding to marked absorption fea- 
tures common to all lines. These broad intervals (typically 5 to 
20 km s ') are given in the two first columns of Table[3] The col- 

* The uncertainties on the moments of the distributions are the statis- 
tical standard errors computed for the first- and second-order moments 
as (t/ and cr/ V2A', respectively, where cr is the standard deviation 
and N the size of the sample. 
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umn densities given in columns 4, 5, and 6 of Table[3]are derived 
assuming a single excitation temperature of 4.4 K for '^CH^ and 
'^CH+, and of 3.0 K for SH+ (see AppendixH). Following the 
method set in Sect. 13.21 to keep only the most reliable absorp- 
tion features, a 3cr detection level was adopted. Throughout, ev- 
ery measurement below this level is considered as a lower limit. 
Conversely, we adopt a conservative lower limit o f 2.3 on the 
optical depth for the saturated velocity intervals tNeufeld et alj 
I2010bh . 

4.1. Variation oftlie '^C/'^C isotopic ratio across tlie Galactic 
disk 
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ratio and A^('^CH^) computed from Eq. |4] for the veloc- 
ity intervals towards DR21(OH), G34.3-H0.1, W31C, W33A, 
W49N, and W51 are given in columns 9 and 10 of Table [3] 
Towards SgrA*+50 and SgrB2(N), because most of the gas 
appe ars to be associated with the Galactic centre environment 
(Ro driguez-Fernandez et al.l l2006h . a '^C/'^C abundance ratio 
of 20 is assumed everywhere except for the velocity interval -20 
to -1-30 km s a velocity domain where the absorption features 
are also associated with gas in the Galactic plane, and where we 
use two alternative values, 20 and 60, to bracket the result. 

The uncertainties given in Eq. |4] correspond to the standard 
deviation of the best least-squares fit performed by Milam et al. 
(2005) on the CN, CO, and H2CO data. They do not take into 
account the uncertainty on Re due to the random motion of inter- 
stellar clouds. To do so, we relied on the analysis of the HI emis- 
sion in the first Galactic quadrant by Elmegreen & Elmegre^ 
(1987): they found that a large part of the mass of the diffuse gas 
is distributed into ~ 200 pc superclouds, separated along spiral 
arms by 1.5 kpc. Their one-dimensional internal velocity dis- 
persion of 5.3 km s ' is the main source of uncertainty on Re- 
Taking into account the (/, v) location of the absorbing gas ob- 
served towards each source (see the equation of Rq in the caption 
of Table O, we obtain a maximal uncertainty of about 30 % on 
Re- When combined with the errors given in Eq. |4] we find a 
total uncertainty of about 50% on the '^C/'-'C ratio and the sub- 
sequent A^('^CH+) given in columns 9 and 10 of Table[3] 

4.2. Comparison of the column densities of CH+ and SH"^ 
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Figure 6. N(^^C)/NC^C) column density ratio as a function 
of the galactocentric distance Rq. The red squares are from 
this work, and therefore evaluated with measurements of the 
NC^CH^)/NC^CH^) column density ratio. The blue circles 
are from previous measurements of the A^('^CN)/A^('^CN), 
NC^C 0)/NC^C0) and MH^^CO)/MHf CO) column density 
ratios (*Milam et al.ll2005l and references therein) while the green 
curve corresponds to a Unear least-squares fit applied to these 
data. 



In Fig.|6] we display the '^C/'^C column density ratio de- 
rived from the present '-CH^ and '^CH^ data as a function of the 
Galactocentric distance Rq (given in columns 7 and 8 of Table 
O which was computed assuming a flat Galactic rotation curve. 
We compared these results to those deduced from previous ob- 
servations of CN, CO, H2CO and their respective isotopologues 
(•Milam et al. 2005). We found that the five firm values (indicated 
in Table 3) derived from the simultaneous detections of '^CH^ 
and '^CH^ absorption lines over the same velocity range are con- 
sistent with those derived from the neutral species. The 25 lower 
limits inferred from sa turated '^CH'^ Une s (dashed symbols) are 
also consistent with Milam et"aD (l2005h . This result not only 
suggests that the isotopic ratios measured with ions and neu- 
trals are not substantially influenced by chemical fractionation 
processes, but it also validates the use of the empirical relation 



\2r^ /I3 



C/'^C = 6.2(+1.0)/?G + 18.7(+7.4) 



(4) 



found by iMilam et all (|2005|) to infer the ^^CH^ column den- 
sities from those measured in the '-^CH^ spectra. The '^C/'"^C 
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Figure?. SH^ column density as a function of the CH^ col- 
umn density per broad velocity interval (see Table |3]l. The open 
squares and open cir cles are from the present analysis and from 
iMenten et al.l (1201 1), respectively. The blue and purple points 
are from the absorption lines observed along the Galactic cen- 
ti-e sight lines (SgrA*4-50, SgrB2(N), and SgrB2(M)) and along 
the other sight lines (DR21(OH), G34.3H-0.1, W31C, W33A, 
W49N, and W51), respectively. The black dashed, dotted, and 
dashed-dotted curves indicate A^(SH+)/A^(CH+) ratios of 0.01, 
0.1, and 1, respectively. 



Fig. |2]displays the SH^ and CH^ column densities inferred 
for each broad velocity interval. The CH^ column densities 
were derived either from the '^CH+ profile, where unsaturated. 
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or from '^CH^ and the above isotopic ratio in the other case. 
The data set includes the results of the present study and those 
obtained b y Menten et al. (2011) towards SgrB2(M). While 
iMenten et alJ (1201 ll) reported the detection of two absorption 
components at A?(CH+) ~ 3x10'^ andA^(SH+) ~ 10'^ cm^^, we 
removed those points from Fig. Q because the ^^CH"^ spectrum 
in the corresponding velocity range clearly exhibits contamina- 
tion of the absorption features by a strong and broad emission 
line. Last, we note that the few points at A^(CH^) < 10'^ cm"^ 
correspond to the faintest '^CH^ absorption features. Since they 
are not detected in SH^, these points weakly constrain our sub- 
sequent analysis of the SHVCH^ratio. 

Taking into account the upper and lower limits shown in Fig. 
|2l we find that A^(CH^) and A^(SH^) span more than two orders 
of magnitude and that the SH^/CH^ column density ratio varies 
from less than 0.01 to 1. Interestingly, the A^(SH+)/A^(CH+) 
ratios observed towards SgrB2(M), SgrB2(N), and SgrA*+50 
are very similar, with a mean value of 0.28 + 0.02, ~ 9 times 
higher than that towards DR21(OH), G34.3H-0.1, W31C, W33A, 
W49N, and W51, <A^(SH+)/A^(CH+)) ~ 0.03 + 0.007. The above 
difference in the mean SH^/CH^ ratio is much larger than the 50 
% uncertainty on the computed '^C/'-'C ratio. Moreover, we ob- 
tain the same difference of values if we use only the few points 
detected in ^^CH^. It is therefore unlikely that this difference can 
be ascribed to an uncertainty on the '^CH^/'-'CH^ ratio. Last, 
jDaflon & Cunha (2004) found that carbon and sulfur have sim- 
ilar abundance gradients across the Galactic disk (with slopes 
of -0.037 and -0.040 dex kpc"', respectively). Consequently, the 
difference observed in the SH^/CH^ column density ratios mea- 
sured on the I ~ and / sight lines is most likely tracing 
variations of both the physical and chemical conditions of the 
diffuse gas sampled in each case. 

Finally, with a correlation coefficient of 0.1, no evident 
chemical relationship seems to stand out from Fig.|7J a surpris- 
ing finding in view of the fact that CH^ and SH^ are clearly 
linked by their dynamics (see Sect. 13.31 ). However, this lack of 
correlation applies to the column densities. In the following, we 
discuss the properties of the abundances relative to hydrogen, a 
discussion that requires the knowledge of A^h- 

4.3. Estimation of the hydrogen coiumn densities in the 
broad veiocity intervals 

To (1) estimate the mean molecular abundances that are to be 
compared with the chemical model predictions (Godard et al., 
in preparation), and (2) to establish a possible relation between 
the CH^ and SH^ abundances, it is essential to estimate the hy- 
drogen column density A^h in the broad velocity intervals over 
which A^(CH^) and A^(SH^) are measured. Since the molecular 
fraction of the gas where CH^ and SH^ are detected is low (0.4 
on average, see AppendixO. both atomic and molecular hydro- 
gen HI and H2 are needed to estimate A^h- 

The method consists in separately evaluating the column 
densities of HI and H2, using t he VLA observations of the 

[al] 
tail 



A21 cm absorption Une of H (iKod 119971: iFish et 



iDwarakanath et alJl2004l: iPandian et alJl2008t iLang et 



and a relevant tracer for the molecular hydrogen because H2 
is not directly observable. Then, A^h = N(HI) + 2A^(H2). In 
Appendix ICl we discuss the validity of using the HIFI observa- 
tions of CH and HF to compute the H2 column densities. If avail- 
able, HF is preferentially used in the following section to infer 
A^(H2) and the ensuing CH+, '^CH+,andSH+ mean abundances. 
If not, A^(H2) is derived from CH, assuming a HF/CH mean 
abundance ratio of 0.4 (a value defined with a large standard de- 



viation of 0.25), deduced from columns 5 and 6 of Table|4] Last, 
we compare these values of A^h to those inferred, as in Godard et 
al. (20 10), from the an alysis of the 2MASS survey (Cutri et al] 
I2003h . iMarshall et al.l (l2006h have measured the near infrared 
colour excess in large areas of the inner Galaxy (|/| < 100°, 
\b\ < 10°) to obtain the visible extinctions (Ay ~ lOA^:), pro- 
viding an estimate of the total hydrogen column density along 
the lines of sight. However, because of the low resolution of the 
2MASS extinction analysis (~ 15 arcmin), the uncertainty on 
A^H (computed as the standard deviation of the extinction mea- 
sured along the four closest lines of sight surrounding a given 
source) can be important, and is as high as 50% for DR21(OH), 
W31C, and W33A. The last columns of Table|4]show that these 
two independent measurements of A^h differ by less than 25 %. 



5. Results and discussion 

5. 1 . Chemicai properties of the gas seen in absorption 
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Figure 8. A^(SH^) /A^(CH^) column density ratio as a function of 
the CH^ mean abundance. The CH+, SH^, and H column den- 
sities are computed from CH+, i^CH+, SH+, HI, HF, and CH 
opacities integrated over the velocity intervals given in Tables 
|3] and |4] The blue and purple points are from the absorption 
lines observed along the Galactic centre sight lines (SgrA*+50 
and SgrB2(N)) and along the other sight hnes (DR21(OH), 
G34.3H-0.1, W31C, W33A, W49N, and W51), respectively. The 
red line corresponds to a least-squares fit of the latter data. 



In Fig.m we display the A^(SH+)/A?(CH+) column density 
ratio as a function of the CH^ mean abundance (with respect to 
the total hydrogen column density A^h) in tegrated over t he ve- 
locity intervals given in Table |4] Because 'Me nten et al.l (1201 ih 
used a different method than we did to estimate the H2 column 
densities Q the data obtained towards SgrB2(M) are not included 
in this plot. We found that both the mean abundances and the 
abundance ratios vary by two orders of magnitude in the diffuse 
ISM sampled by those lines of sight. In addition. Fig. |8]reveals 
major differences between the results obtained on the / ~ and 



' In lMenten et"aD j201ll) the H2 column densities are deduced from 
those of HCO+, assuming Af(HC0+)/A^(H2) = 5 x 10 '' dLucas & Liszj 
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the I + sight lines. While the CH"^ mean abundances span the 
same range of values in both cases, the SH^/CH^ ratio mea- 
sured towards SgrA*+50 and SgrB2(N) shows no correlation 
with A^(CH+)/A^H- Conversely, the points corresponding to the 
observations performed on the / + sight lines exhibit a trend 



A^(SH+)/A?(CH+) ~ 0.09 [lO^ x A?(CH+)/A?h] 



-1.4 



(5) 



with a correlation coefficient of 0.8. These are the results that 
have to be compared with the predictions of chemical models 
applied to the diffuse ISM. 

5.2. Carbon and sulfur chemistries in the diffuse ISM 
5.2.1. UV-driven chemistry 
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Figure 9. Predictions of two PDR models for gas densities rtu - 
50 and 300 cm ^ The CH+ (in red) and SH+ (in blue) relative 
abundances are displayed as functions of the shielding Ay from 
the interstellar radiation field for a slab of gas illuminated on one 
side only. 



In a chemistry entirely driven by the UV-radiation field and 
the cosmic ray particles, the hydrogenation chains of carbon and 
sulfur, and the subsequent productions of CH^ and SH^ are ini- 
tiated by the radiative associations of and with molecular 
hydrogen. 



Hz 
Hz 



CHJ 
SHt 



y and 



(6) 



two reactions with long timescales: 2 x 10'' yr f^^ (SQ cm ^/«h) 

and 1 X 10^yr/- ^(50 cm-V«H), respectively (iHerbslI [1981 
iHerbst et al.|[l989l) . where /h, is the molecular fraction defined 
as = 2«(Hz)/nH- In comparison, CH^ and SH^ are mainly 
destroyed by hydrogenation and dissociative recombinatiorfl 



CH^ + nz 
SH+ + e 



Hz ^ CH+ + H and 



+ H, 



Because tlie hydrogenation of SH+ is liighly endo-energetic: 
AE/k = 6380 K. 



two processes with short timescales: 1 yr f^^^ (50 cm ^ /riuj and 
Syr (r/100K)°-^2(i 33 ^ 10-''/-^.-)(50 cm-'VnH), respectively. 
Because of the lack of efficient production pathways to balance 
their rapid destruction, the CH^ and SH+ abundances predicted 
by UV-dominated chemistry are very low. Fig. |9] displays the 
CH^ and SH+ relative abundances computed with two mod- 
els of PhotoDissociation Regions (PDRjJ illuminated on one 
side as functions of the shielding from the ISRF Ay- For the 
physical conditions of the diffuse gas, we obtain A^(CH^)/A^h ~ 
1.8 X 10"" and A^(SH+)/A^h ~ 1.6 x 10"'\ two to four orders 
of magnitude lower than the observed values, and with a corre- 
sponding abundance ratio never higher than 0.01. 

5.2.2. Alternative models 

In the diffuse interstellar medium sampled by the / sight 
lines, the only production pathways efficient enough to balance 
the fast destruction of CH^ and SH^ are 



C++ Hz 
S+ + Hz 



. CH++H A£/;t = 4640K, 
SH+ + H AE/k = 9860K. 



and 



(8) 
(9) 



Since these reactions are highly endothermic, it has been pro- 
posed that high CH+ and SH+ abundances are the signatures 
of s hock waves p ropagating through the ISM (Draine 1986^ 
Mill ar etal .11 19861) . Comparing the predictions of HD and MHD 
shocks with the C H+ column density obse r ved to wards ^ Oph, 
Draine (1986) and iPineau des Forets et al.l (Il986l) favoured the 
MHD case in which CH+ and SH+ form mainly through 
ion-neutral friction. Including the sulfur chemistry Millar et al] 
(119861) and lPineau des Forets et aLl(ll986h predicted a SH+/CH+ 
abundance ratio increasing from 0.01 to 0.4 for a shock speed in- 
creasing from 9 to 16 km s"', a transverse magnetic field of 5fiG 
and a preshock density of 20 cm""*. These results agree excel- 
lently with our observations. 

Another scenario, alternative to the shock waves, is the TDR 
(turbulent dissipation regions) model. In this model the turbulent 
energy is dissipated in many small-scale magnetized vortices 
in which the ionized and neutral fluids decouple. Dissipation is 
caused by both ion-neutral friction and viscous dissipation at the 
edge of the vortices. Comparing the predictions of TDRs to ob- 
servations of CH, CH+, OH, and HCO+ in the local diffuse gas, 
Godard et al. (2009) also favoured models in which the dissipa- 
tion is dominated by the ion-neutral friction and where the pro- 
duction of CH+ and SH+ via reactions [8] and |9] is triggered by 
the ambipolar diffusion. An analysis of the SH+/CH+ ratio ob- 
tained in the framework of the TDR model will be presented in 
a forthcoming paper (Godard et al. in prep.). 

While the scenarios of shocks and vortices could equally ap- 
ply to gas sampled by the I ~ sight lines, an alternative chem- 
ical process may be at work there. Large amounts of the diffuse 
gas detected along SgrA*+50, SgrB2(N), and SgrB2(M) belong 
to the Central Molecular Zone (CMZ). Because the CMZ is per- 
vaded by a strong X-ray radiation field, high CH+ and SH+ abun- 
dances could be due to regions where C++ and S++ ions co-exist 
with Hz and form CH+ and SH+ by the reactions 



■ H? 



CH+ + H+ 



and 



(7) S++ + Hz SH+ + H+ 



(10) 
(11) 



' The Meudon PDR model employs a one-dimensional chemical 
code in which a slab of gas with a given density profile is ill uminated 
by the ambient interstellar radiation field jLe Petit et all 200^ . 
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as proposed by iLanged ( 19781). U s ing the high rate of reaction 
(HB obtained bv lChen et al.l(l2003h . lAbel et al.l (12008 !) found that 



the predicted SH* column density for molecular gas surrounding 
an active galactic nucleus is two orders of magnitude higher than 
that predicted by UV-dominated chemistry. Since rate measure- 
ments of reaction (fTOb point to lower values, this chemical pro- 
cess could account for the high mean SH^/CH^ abundance ra- 
tio derived towards SgrA*-i-50, SgrB2(N), and SgrB2(M) (about 
nine times higher than that obtained along the other Galactic 
sight lines, see Sect. l4.2l i. 

The above results are in line with the detection of large 
amounts of warm diffuse gas recently identified with the 3fim 
absorption lines of H+ (Oka et.al. 2005; Geballe & O ka 2010) in 
the CMZ. The temperature (~ 200 - 300 K) of the low-density 
gas phase is found to be considerably higher than t hat of typi- 
cal di ffuse clouds, and unique to the Galactic CMZ dGoto et al.l 
120081) . They also coiToborate the specific dynamics of molec- 
ular clouds associated with the Galactic centre regions (e.g. 
iRodriguez-Fernandez et al ]l2006h . 

6. Summary and perspectives 

We have presented the analysis of HerschelfHlFl observations 
of the ground-state transitions of CH^, '^CH^, and SH^, all de- 
tected in absorption against the submillimetre dust continuum 
of distant star-forming regions and the Galactic centre sources, 
SgrA*-i-50 and SgrB2(N). The velocity range over which the ab- 
sorption features are detected corresponds to diffuse or translus- 
cent environments. The deconvolution of the hyperfine structure 
embedded in the SH^ 1,2 - 0, 1 spectra, and the independent 
decomposition of the absorption domains in Gaussian velocity 
components allowed us to identify many velocity components 
per sight line and to perform a cross comparison of the dynami- 
cal and chemical signatures of those three species. 

This study provides the following main results. (1) The 
linewidth distributions of CH^, '^CH^, and SH^ are found to 
be similar and likely trace the kinematics of the chemical pro- 
duction processes of these species convolved with that of the 
turbulent and Galactic dynamics of the diffuse ISM. (2) These 
lines are broad (~ 4.2 km s"'), similar to those found in visible 
absorption lines in the solar neighbourhood, and broader than 
those of HCO^ and CN along the same lines of sight. (3) The 
SH^/CH^ abundance ratio covers a broad range of values from 
0.01 to more than 1, shows higher values in warmer environ- 
ments (such as the Galactic centre clouds), and appears to be 
proportional to (A^(CH^)/A^h) ' "^ in the diffuse gas sampled by 
the I + sight lines. (4) As for CH^, the SH^ abundances cannot 
be reproduced by UV-driven chemistry in the diffuse gas. 

The unique properties of the carbon and sulfur chemistries 
support the framework of a warm chemistry triggered by turbu- 
lent dissipation (either in shocks or intense velocity shears) that 
selectively enhances the production of SH^ for which the forma- 
tion endothermicity is the highest. A detailed comparison of the 
TDR model predictions with these observational results will be 
given in a forthcoming paper (Godard et al. in prep.). 

Acknowledgements. We are most grateful to the referee for providing construc- 
tive comments and helping in improving the content of this paper HIFI has been 
designed and built by a consortium of institutes and university departments from 
across Europe, Canada and the United States (NASA) under the leadership of 
SRON, Netherlands Institute for Space Research, Groningen, The Netherlands, 
and with major contributions from Germany, France and the US. Consortium 
members ai'e : Canada: CSA, U. Waterloo; France : CESR, LAB, LERMA, 
IRAM; Germany : KOSMA, MPIfR, MPS; Ireland : NUI Maynooth; Italy : 
ASI, IFSI-INAF, Osservatorio Astrofisico di Ai-cetri-INAF; Netherlands : SRON, 
TUD; Poland : CAMK, CBK; Spain : Observatorio Astronomico Nacional 



(IGN), Centro de Astrobiologia; Sweden : Chalmers University of Technology 
- MC2, RSS & GARD, Onsala Space Observatory, Swedish National Space 
Boai'd, Stockholm University - Stockholm Observatory; Switzerland : ETH 
Zurich. FHNW; USA : CalTech, JPL, NHSC. BG, EE, MG, and MDL acknowl- 
edge the support from the Centre National de Recherche Spatiale (CNES), and 
from ANR through the SCHISM project (ANR-09-BLAN-231). BG, JC, and 
JRG thank the Spanish MICINN for funding support through grants, AYA2009- 
07304 and CSD2009-00038 



References 

Abel, N. R, Federman, S. R., & Standi, R C. 2008. ApJ, 675, L81 
Abia, C, Cunha, K., Cristallo, S., et al. 2010, ApJ, 715, L94 
Allen, M. M. 1994, ApJ, 424, 754 
Amano, T. 2010, ApJ, 716, LI 

Brown, A. & Balint-Kurti, G. G. 2000, J. Chem. Phys., 113, 1870 
Chen, D., Gao, H., & Kwong, V. H. 2003, Phys. Rev. A, 68, 052703 
Comito, C. & Schilke, P 2002, A&A, 395, 357 
Crane, P, Lambert, D. L., & Sheffer, Y. 1995, ApJS, 99, 107 
Crawford, I. A. 1995, MNRAS, 277, 458 
Crawford. I. A. & WilHams, D. A. 1997, MNRAS, 291, L53 
Cutri, R. M., Skrutskie, M. F, van Dyk, S., et al. 2003, 2MASS All Sky Catalog 
of point sources., ed. Cutri, R. M., Skrutskie, M. F., van Dyk, S., Beichman, 
C. A., Carpenter J- M.. Chester, T., Cambresy, L., Evans, T., Fowler, J.. Gizis, 
J., Howard, E., Huchra, J., Jarrett, T, Kopan, E. L., Kirkpatrick, J. D., Light, 
R. M., Marsh, K. A., McCallon, H., Schneider, S., Stiening, R., Sykes, M., 
Weinberg, M., Wheaton, W. A., Wheelock, S., & Zacarias, N. 
Daflon. S. & Cunha, K. 2004. ApJ, 617, 1115 

Danks, A. C, Federman, S. R., & Lambert, D. L. 1984, A&A, 130, 62 

Dowell, C. D., Lis, D. C, Serabyn, E., et al. 1999, in Astronomical Society of 
the Pacific Conference Series, Vol. 186, The Central Parsecs of the Galaxy, 
ed. H. Falcke, A. Cotera, W. J. Duschl, F. Melia, & M. J. Rieke, 453 

Draine, B. T. 1986, ApJ, 310, 408 

Draine, B. T. & Katz, N. 1986. ApJ, 310, 392 

Dwarakanath, K. S., Goss, W. M., Zhao, J. H., & Lang, C. C. 2004, Journal of 

Astrophysics and Astronomy, 25, 129 
Elmegreen, B. G. & Elmegreen, D. M. 1987, ApJ, 320, 182 
Falgarone, E., Godard, B., Cernicharo, J., et al. 2010, A&A, 521, L15+ 
Falgarone, E., Pineau des Forets, G., & Roueff, E. 1995, A&A, 300, 870 
Federman, S. R., Rawlings, J. M. C, Taylor, S. D., & Williams, D. A. 1996, 

MNRAS, 279, L41 
Federman, S. R., Strom, C. J., Lambert, D. L., et al. 1994, ApJ, 424, 772 
Fish, V. L., Reid, M. J., Wilner, D. J., & Churchwell, E. 2003, ApJ, 587, 701 
Geballe, T. R. & Oka, T. 2010, ApJ, 709, L70 
Gerin, M., de Luca, M., Black, J., et al. 2010a, A&A, 518, L110+ 
Gerin, M., de Luca, M., Goicoechea, J. R., et al. 2010b, A&A, 521, L16+ 
Godard, B., Falgaione, E., Gerin, M., Hily-Blant, P, & de Luca, M. 2010, A&A, 

520, A20+ 

Godard, B., Falgarone, E., & Pineau Des Forets, G. 2009, A&A, 495, 847 
Goto, M., Usuda. T., Nagata, T., et al. 2008, ApJ, 688, 306 
Greaves, J. S. & WilUams, P G. 1994, A&A, 290, 259 
Gredel, R. 1997, A&A, 320, 929 

Gredel, R., van Dishoeck, E. F, & Black, J. H. 1993, A&A, 269, 477 
Gry, C, Boulanger, F, Nehme, C, et al. 2002, A&A, 391, 675 
Hammami. K., Owono Owono, L. C, & Stauber. R 2009, A&A. 507, 1083 
Haud, U. & Kalberla. P M. W. 2007, A&A, 466, 555 
Herbst, E. 1985, ApJ, 291, 226 

Herbst, E., Defrees, D. J., & Koch, W. 1989, MNRAS, 237, 1057 
Indriolo, N., Oka, T., Geballe, T. R.. & McCall, B. J. 2010, ApJ, 711. 1338 
Jenniskens, P, Ehi-enfreund, R, & Desert, F 1992, A&A, 265, LI 
Joulain, K., Falgarone, E., Pineau des Forets, G., & Flower, D. 1998, A&A, 340, 
241 

Kerr, F J. & Lynden-Bell, D. 1986, MNRAS, 221, 1023 
Koo, B. 1997, ApJS, 108, 489 

Lacour, S., Ziskin, V., Hebrard, G., et al. 2005, ApJ, 627, 251 

Lambert, D. L., Sheffer, ¥., & Crane, R 1990, ApJ, 359, L19 

Lang, C. C, Goss, W. M., Cyganowski, C, & Clubb, K. I. 2010, ApJS, 191, 275 

Langer, W. D. 1978, ApJ, 225, 860 

Le Petit, F, Nehme, C, Le Bourlot, J., & Roueff, E. 2006, ApJS, 164, 506 

Lesaffre, P, Gerin, M., & Hennebelle, R 2007, A&A, 469, 949 

Lim, A. J.. Rabadan, I., & Tennyson, J. 1999, MNRAS, 306, 473 

Liszt, H. & Lucas, R. 2002, A&A, 391, 693 

Liszt, H. S., Lucas, R., & Pety, J. 2006, A&A, 448, 253 

Lucas, R. & Liszt, H. 1996, A&A, 307, 237 

Magnani, L. & Salzer, J. J. 1989, AJ, 98, 926 

Magnani, L. & Salzer, J. J. 1991, AJ, 101, 1429 



9 



B. Godard et al.: Comparative stiidy of CH* and SH* absorption lines. 

Maier, J. P., LaMn, N. M., Walker, G. A. H., & Bohlender, D. A. 2001, ApJ, 553, 

267 

Marshall, D. J., Robin, A. C, Reyl6, C, Schultheis, M., & Picaud, S. 2006, 

A&A, 453, 635 
Mattila, K. 1986, A&A, 160, 157 

Menten, K. M., Wyrowski, F., Belloche, A., et al. 201 1, A&A, 525, A77+ 
Milam, S. N., Savage, C, Brewster, M. A., Ziurys, L. M., & WyckoflF, S. 2005, 
ApJ, 634, 1126 

Millar, T. J., Adams, N. G., Smith, D., Lindinger, W., & ViUinger, H. 1986, 

MNRAS, 221,673 
Millar, T. J. & Hobbs, L. M. 1988, MNRAS, 231, 953 

Neufeld, D. A., Goicoechea, J. R., Sonnentrucker, R, et al. 2010a, A&A, 521, 
L10+ 

Neufeld, D. A., Kaufman, M. J., Goldsmith, R R, Hollenbach, D. J., & Plume, 

R. 2002, ApJ, 580. 278 
Neufeld, D. A., Sonnentrucker, R, Phillips, T. G., et al. 2010b, A&A, 518, L108+ 
Neufeld, D. A. & Wolfire, M. G. 2009, ApJ, 706, 1594 
Neufeld, D. A., Wolfire, M. G., & Schilke, R 2005, ApJ, 628, 260 
Novotny, O., Mitchell, J. B. A., LeGarrec, J. L., et al. 2005, Journal of Physics 

B Atomic Molecular Physics, 38, 1471 
Nyman, L. & Millar, T. J. 1989, A&A, 222, 231 

Oka, T., Geballe, T. R., Goto, M., Usuda, T., & McCall, B. J. 2005, ApJ, 632, 
882 

Ott, S., Science Centre, H., & Space Agency, E. 2010, ArXiv e-prints 

Pan, K., Federman, S. R., Cunha, K., Smith, V. V., & Welty, D. E. 2004, ApJS, 

151, 313 

Pan, K., Federman, S. R., Shefter, Y, & Andersson, B. 2005, ApJ, 633, 986 
Pandian, J. D., Momjian, E., & Goldsmith, P R 2008, A&A, 486, 191 
Persson, C. M., Black, J. H., Cemicharo, J., et al. 2010, A&A, 521, L45+ 
Pety, J. 2005, in SF2A-2005: Semaine de I'Astrophysique Francaise, ed. 

F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 721 — I- 
Pineau des Forets, G., Roueft', E., & Flower, D. R. 1986, MNRAS, 223, 743 
Plume, R., Kaufman, M. J., Neufeld, D. A., et al. 2004, ApJ, 605, 247 
Rachford, B. L., Snow, T. R, Tumlinson, J., et al. 2002, ApJ, 577, 221 
Reach, W. T., Dwek, E., Fixsen, D. J., et al. 1995, ApJ, 451, 188 
Rodriguez-Fernandez, N. J., Combes, F., Martin-Pintado, J., Wilson, T. L., & 

Apponi, A. 2006, A&A, 455, 963 
Roelfsema, R R., Helmich, F R, Teyssier, D., et al. 2012, A&A, 537, A17 
Rudolph, A. L., Fich, M., Bell, G. R., et al. 2006, ApJS, 162, 346 
Savage, B. D., Bohlin, R. C, Drake, J. F, & Budich, W. 1977, ApJ, 216, 291 
Savage, B. D. & Sembach, K. R. 1996, ARA&A, 34, 279 
Savage, C, Apponi, A. J., & Ziurys, L. M. 2004, ApJ, 608, L73 
Shefl'er, Y., Rogers, M., Federman, S. R., et al. 2008, ApJ, 687, 1075 
Shull, J. M. & Beckwith, S. 1982, ARA&A, 20, 163 
Snow, T. R, Destree, J. D., & Jensen, A. G. 2007, ApJ, 655, 285 
Snow, T. R & McCall, B. J. 2006, ARA&A, 44, 367 
Snow, T. R, Rachford, B. L., Tumlinson, J., et al. 2000, ApJ, 538, L65 
Sofia, U. J. & Meyer, D. M. 2001, ApJ, 554, L221 

Sonnentrucker, P, Neufeld, D. A., Phillips, T. G., et al. 2010, A&A, 521, L12-H 
Thorbum, J. A., Hobbs, L. M., McCall, B. J., et al. 2003, ApJ, 584, 339 
Weselak, T., Galazutdinov, G. A., Musaev, F. A., & Krelowski, J. 2008, A&A, 

484, 381 

Xie, T., Allen, M., & Langer, W. D. 1995, ApJ, 440, 674 

Zhu, C, Krems, R., Dalgamo, A., & Balakrishnan, N. 2002, ApJ, 577, 795 



10 



B. Godard et al.: Comparative study of CH* and SH* absorption lines. 




-300 -200 -100 100 200 300 



-200 -100 100 200 300 400 



3.00 




240 



240 



180 60 

LSR velocity (km s"^) 



240 



Figure 2. Spectra observed at the frequencies of the ground-state transitions of CH^ (black), '^CH^ (red), and SH^ (blue) in the 
direction of DR21(OH), G34.3-H0.1, W31C, W33A, W49N, W51, SgrA*+50, and SgrB2(N). All spectra have been normaUzed to 
the continuum temperatures. To make the absorption features discernible, the '^CH^ and SH^ signals were mutiplied by a factor of 
2 towards SgrA*-i-50, and SgrB2(N), and by factors of 2 and 3, respectively, towards the other sources. For more clarity the '^CH^ 
and SH^ spectra were shifted from the CH^ spectra by 0.5 and 1.0. 
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Figures. Absorption profiles observed in the direction of DR21(OH), G34.3-H0.1, W31C, W33A, W49N, W51, SgrA*-H50, and 
SgrB2(N) in the ground-state transitions of CH^ (black) and '^CH^ (red), and in the main hyperfine component (see Tab.|2]i of the 
ground-state transition of SH^ (blue). All spectra have been normalized to the continuum temperatures. To make the absorption 
features discernible, the '^CH^ and SH^ signals were mutiplied by a factor of 2 towards SgrA*-i-50, and SgrB2(N), and by factors 
of 2 and 3, respectively, towards the other sources. For more clarity the '^CH^ and SH^ spectra were shifted from the CH^ spectra 
by 0.5 and 1.0. The x- and y- scales were chosen to display the velocity structure of the absorption features in detail. 
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Figure 4. Result of the hyperfine structure deconvolution code applied on the SH^ (1,2 <— 0, 1) absorption spectrum observed 
towards SgrB2(N) and W49N. The top panel displays the original data and the resulting decomposition in three hyperfine com- 
ponents, aligned to the same velocity scale. The hyperfine components are shifted by their velocity shift Ay/, (see Table [2]i in the 
bottom panels, to show how they can be combined (green line) to reproduce the original data (black Une). For more clarity, the 
hyperfine components were vertically shifted from the original data. 
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Figures. Top: (a) and (b) histogram of the Hnewidths of the CH+, '^CH+, SH+ (this work), and HCO+ (iGodard et al.ll2010h 
velocity components obtained with the multi-Gaussian decomposition procedure (see Table lA.lb ; (c) linewidth distributions of 
CH^ and CH observed by Crane et al. (1995) in the local diffuse medium. Bottom: first- (mean) and second- (standard devia- 
tion) order moments of the Av distributions issued from the combined CH^, '-^CH^, and SH^ data (red squares) and from the 
CH^ data observed in the solar neighbourhood (blue circles); the moments are computed over increasing Av intervals: mean = 

J^""" AvNiAv)d(Av)/ J^^""'" N(Av)d(Av); standard deviation - [^'^'^''""(Al' - meanfNiAv)diAv)/ J^"""^ N(Av)d(Av)y'^, where 
N{Av) is the linewidth distribution. 
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Table 3. Comparison between the column densities of CH^, '''CH^, and SH^obtained with the integrations of the respective opacities over several 
velocity intervals. The two last columns contain the ratios = A'(SH+)/A'(CH+) and ri^ = jV(S H+)/A^("CH+ ) x '^C /'-C. In column 10, CH+ 
column densities are inferred from those of '3CH+ and the '^c/'^C Galactic gradient derived bv lMilametaP d2005h : A'(CH+) = nC^CH^) x 
'^C/'^C. The blank spaces correspond to unusable '^CH+ or SH+ data owing to the contamination of an emission line from the SFRs, or because 
the spectrum measured over the associated velocity interval could not be deconvolved from the SH^ hyperfine structure. 
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6.2 


6.6 


59 


1.8± 1.1 


< 


3.0 




<0.4 




18 


25 




> 3.8 


1.7± 1.0 


1.5 ±0.2 


5.6 


6.2 


56 


8.1±2.1 


< 


0.4 




0.2 ± 


0.1 


25 


31 




> 4.2 


2.2 ± 0.9 


2.8 ± 0.2 


5.2 


5.6 


52 


14.7 ± 2.7 


< 


0.5 




0.2 ± 


0.1 


31 


35 




> 3.0 


4.7 ± 0.8 


1.6± 0.2 


5.0 


5.2 


50 


8.0± 1.7 


< 


1.6 




0.6 ± 


0.2 


35 


39 


E 


> 3.0 


> 1.2 


> 1.1 


4.7 


5.0 


49 


> 5.3 













* E = absorption line profile observed in the star-forming region; T^,, may be underestimated, hence the lower limit on A'(CH+), A'('^^CH+), and 
A'(SH^). 1,2,3,4,5 = velocity intervals where ''CH^ and "CH^ are observed simultaneously with a 3(T detection level. 

" Galactocentric distances derived as Rr, = Tr-^r-^^^^TT— rrr with Rn = 8.5 kpc and On = 220 km s"' as recommended by the lAU 

^Isr+'^O sin(/)cos(o) u f u J 

jKerr & Lvnden-Belllll986 D. and assuming a flat Galactic rotation curve {d(R(;) = Oq). 

* Towards SgrA*+50 and SgrB2(N), when the comparison is possible, A'(CH+) derived from '^CH+ and '^^CH+ and the ratios ri2 and agree 
within the combination of their respective errors and the uncertainty of 50 % on the '-C/"C ratio (see main text). 
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Table 3. continued. 



^min 










iV(SH+) 


A'(''CH+) 


R" 

e.min 


Da 




A'(CH+)'' 






r^3' 




(km S" 




(kms"') 




(10'3 cm-2) 




(10'2 cm-2) 


(10'^ cm-2) 


(kpc) (kpc) 




(10" cm-2) 


XlO 




XlO 






39 


47 


E 


> 3.4 


> 


5.3 


> 1.9 


4.3 


4.7 


47 


> 8.8 












W49N 




-8 


2 


E 


> 0.5 








8.4 


8.5 


71 
















5 


20 


E 


> 10.4 


> 


19.6 


> 6.6 


7.5 


8.2 


68 


>44.8 














20 


30 


4 


5.1 ±0.1 




2.6 ±0.5 


0.7 ±0.1 


7.1 


7.5 


64 


4.3 ± 1.2 




0.5 ± 


0.1 


0.6 ± 


0.2 




30 


37 




> 5.0 




2.3 ±0.5 


1.3± 0.1 


6.8 


7.1 


62 


8.1 ± 1.5 


< 


0.4 




0.3 ± 


0.1 




37 


44 




> 4.7 




2.6 ± 0.4 


1.1± 0.1 


6.6 


6.8 


60 


6.9 ± 1.4 


< 


0.6 




0.4 ± 


0.1 




44 


49 




> 3.1 


< 


0.4 


0.7 ± 0.1 


6.4 


6.6 


59 


4.2 ± 1.0 


< 


0.1 




<0.1 






49 


54 




> 3.0 




2.0 ± 0.4 


0.6 ± 0.1 


6.3 


6.4 


58 


3.7 ±0.9 


< 


0.7 




0.6 ± 


0.2 




54 


72 




> 11.9 




6.2 ±0.7 


3.6±0.2 


5.7 


6.3 


56 


20.3 ±3.1 


< 


0.5 




0.3 ± 


0.1 




72 


80 




0.4 ±0.02 


< 


0.5 


< 0.1 


5.6 


5.7 


54 


< 0.6 


< 


1.2 








W51 




-1 


3 




0.2 ± 0.01 








8.4 


8.5 


71 
















3 


11 




1.3 ± 0.02 








8.0 


8.4 


69 
















11 


16 




0.3 ±0.01 


< 


0.5 




7.8 


8.0 


68 




< 


1.9 










16 


20 




0.2 ± 0.01 


< 


0.5 




7.6 


7.8 


66 




< 


3.2 










20 


32 




0.6 ±0.02 


< 


0.8 


< 0.2 


7.1 


7.6 


64 


< 1.1 


< 


1.3 










36 


43 




0.3 ±0.01 




1.1 ± 0.6 


0.4± 0.1 


6.8 


7.0 


61 


2.6± 1.2 




3.9 ± 


2.3 


0.4 ± 


0.3 




43 


50 




2.5 ± 0.05 




3.0 ±0.7 


0.4± 0.1 


6.5 


6.8 


60 


2.5± 1.1 




1.2± 


0.3 


1.2± 


0.6 




50 


60 


5 


4.5 ± 0.1 




2.7 ± 0.7 


2.0 ± 0.2 


6.3 


6.5 


59 


11.9± 2.3 




0.6 ± 


0.2 


0.2 ± 


0.1 




60 


69 


E 


> 5.4 


> 


2.1 


> 2.5 


6.0 


6.3 


57 


> 14.4 














69 


75 


E 


> 0.8 






> 0.6 


5.9 


6.0 


56 


> 3.5 













* E = absorption line profile observed in the star-forming region; Tex may be underestimated, hence the lower limit on A'(CH+), A'('^^CH+), and 
A'(SH^). 1,2,3,4,5 = velocity intervals where ''CH^ and '^CH^ are observed simultaneously with a 3(T detection level. 

" Galactocentric distances derived as Rc, = -Kn -ir^^r^^^^^T^r^^rr: with Rn = 8.5 kpc and Qn = 220 km s"' as recommended by the lAU 
jKerr & Lvnden-Be ll 1986), and assuming a flat Galactic rotation curve (6(Rg) = do)- 

* Towards SgrA*+50 and SgrB2(N), when the comparison is possible, A'(CH^) derived from '^CH^ and '^CH^ and the ratios rn and ri^, agree 
within the combination of their respective errors and the uncertainty of 50 % on the '-C/'^C ratio (see main text). 



16 



B. Godard et al.: Comparative study of CH* and SH^ absorption lines. 



Table 4. HI, H2, and total hydrogen column densities per velocity interval. In the last column, the total hydrogen column densities are estimated 
over the entire sight line. 









A'(H)" 


A'(CH)'' 






N(H2Y 












r 

(cm -) 


(kms-') 


(kms-') 




(cm-2) 


(cm-2) 


(cm-2) 




(cm--) 


(cm--) 




(cm-2) 


(cm-2) 










xlO^o 


xlO'"* 


xlO'2 




xlO^« 


xl02« 




xlO^" 


xlO^o 




xl02« 


SgrA*+50 


-168 


-150 




1.5 






















-150 


-123 




4.3 


21.4±3.3 


35.3 + 


1.5 


24.4+ 3.8 


10.1 + 


0.4 


53.2+7.6 


24.5 + 


0.9 




-123 


-95 




1.9 


8.0+ 1.8 


23.5 + 


1.1 


9.2+ 2.1 


6.7 + 


0.3 


20.3 + 4.2 


15.3 + 


0.6 




-95 


-75 




1.9 


4.9 ±0.3 


6.9 + 


0.7 


5.7+0.4 


2.0 + 


0.2 


13.2+0.7 


5.8 + 


0.4 




-75 


-45 




24.7 


17.4 ± 2.0 


> 62.3 




19.8+2.2 


> 17.8 




64.4+4.5 


> 60.3 






-45 


-15 




14.3 


18.9± 1.8 


> 63.9 




21.6+ 2.1 


> 18.3 




57.5 + 4.2 


> 50.8 




480 


-15 







47.7 


22.8 + 4.7 


> 61.2 




26.0+5.4 


> 17.5 




99.7+ 10.8 


> 82.7 









26 




61.8 


11.2+2.3 


> 26.2 




12.8+2.7 


> 7.5 




87.5+5.3 


> 76.8 






26 


49 


E 


>39.0 




> 101.3 






>28.9 






> 96.9 






49 


67 


E 


>31.0 




> 35.0 






>10.0 






> 51.0 






67 


75 


E 


> 8.7 




> 0.9 






> 0.2 






> 9.2 






SgrB2(N) 


-135 


-100 




1.7 




26.0 + 


0.8 




7.4 + 


0.22 




16.6 + 


0.4 




-82 


-62 




5.4 




46.7 + 


1.3 




13.4 + 


0.38 




32.1 + 


0.8 




-62 


-35 




17.4 




> 111.8 






>32.0 






> 81.3 






-35 


-15 




9.2 




82.1 + 


2.8 




23.5 + 


0.80 




56.1 + 


1.6 


450 


-15 


23 




55.8 




138.7 + 


3.6 




39.6 + 


1.01 




135.1 + 


2.0 




23 


38 




1.0 




19.0 + 


0.6 




5.4 + 


0.16 




11.9 + 


0.3 




38 


74 


E 


>69.7 




> 143.0 






>40.9 






> 151.4 






74 


100 


E 


> 14.2 




> 74.6 






>21.3 






> 56.8 






DR21(OH) 


-17 


-9 


E 


>36.3 




> 0.8 






> 0.2 






> 36.8 






-8 


1 


E 


>40.8 


> 12.8 


> 30.9 




> 14.6 


> 8.8 




> 70.1 


> 58.5 






1 


7 




>24.9 


2.0+0.9 


7.0 + 


0.2 


2.3+ 1.0 


2.0 + 


0.05 


> 29.4 


> 28.9 




160 


7 


15 




>38.5 


11.2+0.3 


> 27.4 




12.9+0.3 


> 7.8 




> 64.2 


> 54.1 






18 


26 




>36.3 


< 0.1 


0.3 + 


0.1 


< 0.1 


0.1 + 


0.03 




> 36.5 






G34.3-H0.1 


-4 


5 




5.4 


0.4+0.3 


1.1 + 


0.1 


0.4+ 0.3 


0.3 + 


0.04 


6.2+0.7 


6.0 + 


0.1 




5 


21 




20.1 


4.1 + 0.5 


> 19.8 




4.7+0.6 


> 5.7 




29.4+ 1.2 


> 31.4 






21 


34 




44.2 


4.7+0.5 


> 12.9 




5.3+0.6 


> 3.7 




54.9+ 1.1 


> 51.6 




230 


36 


44 




15.5 


< 0.2 


2.7 + 


0.1 


< 0.2 


0.8 + 


0.04 


< 15.9 


17.1 + 


0.1 




44 


58 


E 


>52.8 


> 9.1 


> 62.9 




> 10.4 


>18.0 




> 73.7 


> 88.7 






58 


66 


E 


>56.1 


> 15.4 


> 37.8 




> 17.6 


> 10.8 




> 91.4 


> 77.7 






W31C 


3 


13 




18.1 


0.9+0.4 


9.4 + 


0.2 


1.0+ 0.4 


2.7 + 


0.1 


20.1 + 0.9 


23.4 + 


0.1 




13 


25 




35.3 


21.8+0.7 


> 59.2 




24.9 + 0.7 


> 16.9 




85.1+ 1.5 


> 69.1 






25 


31 




26.1 


20.1+ 1.3 


> 36.2 




23.0+ 1.4 


> 10.4 




72.1 + 2.9 


> 46.8 






31 


36 




18.5 


9.2+0.4 


> 25.3 




10.5+0.4 


> 7.2 




39.4+0.8 


> 33.0 




220 


36 


44 




36.9 


17.6+0.4 


> 41.0 




20.1 + 0.4 


> 11.7 




77.1 + 0.9 


> 60.3 






44 


51 




1.8 


1.7 + 0.3 


> 8.7 




1.9+0.4 


> 2.5 




5.7+0.7 


> 6.8 






51 


61 




2.5 


0.3 + 0.2 


< 0.2 




0.4+0.3 


< 0.05 




3.3+0.5 


< 2.6 






W33A 


14 


18 




1.7 


0.4+0.1 


0.3 + 0.2 


0.4+ 0.2 


0.1 + 


0.04 


2.5+0.3 


1.9 + 


0.1 




18 


25 




18.0 


2.1 + 0.7 


5.8 + 


0.4 


2.4+ 0.8 


1.7 + 


0.1 


22.9+ 1.5 


21.3 + 


0.2 




25 


31 




22.2 


7.0+0.4 


> 21.1 




8.0+0.4 


> 6.0 




38.3+0.9 


> 34.3 




100 


31 


35 




17.9 


8.1+ 1.0 


> 20.9 




9.2+ 1.2 


> 6.0 




36.3+2.3 


> 29.8 






35 


39 


E 


> 7.7 


> 9.1 


> 20.9 




> 10.4 


> 6.0 




> 28.5 


> 19.6 






39 


47 


E 


> 9.8 


> 5.5 


> 17.9 




> 6.3 


> 5.1 




> 22.4 


> 20.1 






W49N 


5 


20 


E 


>89.1 


> 7.2 


> 75.6 




> 8.2 


>21.6 




> 105.5 


> 132.3 






20 


30 




19.3 


1.6+0.4 


6.1 + 


0.1 


1.8+0.5 


1.7 + 


0.04 


23.0+ 1.0 


22.8 + 


0.1 





* E = absorption line profile observed in the star-forming region ; may be unde r estima t ed, hence the lower limit on A^(HF), and A ^ (CH). 

" From t he observations with the VLA interferometer bv lKool ( Il997h . Ipish et alJ ( l2003h . iPandian et alj j2008h . iDwaralcanath et all l2004h . and 
iLang et al. (2010) a ssumin g a spin temperature of 100 K. 

* From (Gerin et al. 2010b; Gerin et al., priv. comm.). 

From (Neufeld et al. 2010b; Sonnentrucker et al. 2010; Neufeld et al. priv. comm.; Sonnentrucker et al., priv. comm.). 

Estimate from the observation of the CH absorption lines, a mean value of the HF/CH column density ratio of 0.4 (inferred from columns 5 and 
6) and an abundance ratio «(HF)/«(H2) = 3.6 x 10"**, as predicted by UV-dominated chemical models (see Appendix [Ct. 

Estimate from the observation of the HF absorption lines and an abundance ratio ii(HF)/ii(H2) = 3.6 x 10"^, as predicted by UV-dominated 

chemical models (see Appendix ICt. 

^ Estimate of A'h integrated over the entire line of sight from models of the extinction at 2 fim bv lMarshall et al.l j2006h . 
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Table 4. continued. 



V 


u 


* 




A'(CH)* 


A'(HF)'' 


A^(H2)'' 






H , 




H . 




(cm-^) 
xlO^" 


(km s-i) 


(kms"') 




(cm-2) 
xlO^o 


(cm-2) 
xlO" 


(cm-2) 
xlO'2 


(cm-^) 
xlO^o 


(cmr^) 
xlO^o 




(cm -) 
xlO^o 




(cm-2) 
xlO^o 




30 


37 




19.6 


5.2± 1.7 


> 23.7 


5.9± 1.9 


> 6.8 




31.4± 


3.8 


> 33.1 






37 


44 




30.9 


5.8± 1.7 


> 24.4 


6.6 ±2.0 


> 7.0 




44.1 ± 


4.0 


> 44.8 




460 


44 


49 




9.0 


1.0±0.3 


3.3 ±0.1 


1.1± 1.4 


0.9 ± 


0.03 


11.2± 


2.8 


10.8 ± 


0.1 




49 


54 




11.5 


3.2 ± 1.4 


1 1.3 ± 0.2 


3.7± 1.7 


3.2 ± 


0.1 


18.8± 


3.3 


18.0± 


0.1 




54 


72 




52.6 


19.8±3.0 


> 56.0 


22.6 ± 3.4 


>16.0 




97.8 ± 


6.9 


> 84.6 






72 


80 




21.7 


0.8 ±0.4 


0.4 ± 0.1 


0.9 ± 0.4 


0.1 ± 


0.03 


23.5 ± 


0.8 


21.9± 


0.1 




W51 


3 


11 




14.8 


3.1 ± 1.0 


15.2 ±0.4 


3.5± 1.1 


4.3 ± 


0.1 






23.5 ± 


0.2 




11 


16 




4.3 


0.7 ±0.7 


2.1 ±0.1 


0.8 ±0.8 


0.6 ± 


0.03 






5.5 ± 


0.1 




16 


20 




2.3 




0.4 ±0.1 




0.1 ± 


0.03 






2.5 ± 


0.0 




20 


32 




8.6 




1.4±0.1 




0.4 ± 


0.04 






9.4 ± 


0.1 




36 


43 




2.9 




0.1 ±0.1 




0.04 ± 


0.04 






3.0 ± 


0.1 


250 


43 


50 




18.7 


2.3 ±0.7 


> 15.5 


2.6 ±0.8 


> 4.4 








> 27.6 






50 


60 




50.9 




> 41.5 




> 11.8 








> 74.6 






60 


69 


E 


>64.2 


> 12.5 


> 53.1 


> 14.3 


> 15.2 




> 92.8 




> 94.5 






69 


75 


E 


> 10.4 


> 0.8 


> 12.0 


> 0.9 


> 3.4 




> 12.2 




> 17.2 







* E = absorption line profile observed in the star-forming region; Tex may be unde r estima t ed, hence the lower limit on A^(HF), and A ^ (CH). 

° From the observations with the VLA interferometer bv lKod ( Il997h . Ipish et alJ ( l2003h . iPandian et alj |2003) . iDwarakanath et al] j2004h . and 
iLang et al. (2010) assuming a spin temperature of 100 K. 

* From (Gerin et al. 2010b; G erin et al., priv. comm.). 

From (Neufeld et al. 201 Ob ;'Sonnen trucker et al.ll2010l; Neufeld et al. priv. comm.; Sonnentrucker et al., priv. comm.). 

Estimate from the observation of the CH absorption lines, a mean value of the HF/CH column density ratio of 0.4 (inferred from columns 5 and 
6) and an abundance ratio ;!(HF)/;!(H2) = 3.6 x 10"**, as predicted by UV-dominated chemical models (see Appendix [Ct. 

Estimate from the observation of the HE absorption lines and an abundance ratio ii(HF)/ii(H2) = 3.6 x 10 as predicted by UV-dominated 

chemical models (see Appendix ICt. 

^ Estimate of integrated over the entire line of sight from models of the extinction at 2 yjm bv lMarshall etai] ( |200^ . 
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Appendix A: Gaussian decomposition and 
calculation of column densities 



Source: G34 



Source: DR21(0H) 



1.15 
1.10 
1.05 
1.00 
0.95 
0.90 
0.85 

0.00 
-0.10 

1.10 
1.05 
1.00 
0.95 
0.90 
0.85 
0.80 

0.00 
-0.05 





1 : : 




^ SH* (1,2,5/2-0,1,3/2) — ' 







0.00 i^^^^^^^^^^^^^^V^^^^^^^^^^A^ 
-0.05 



-20 



-10 



10 
LSR velocity (km.s"^) 



20 



30 



Figure A.l. Observational data (black lines) compared to 
the multi-Gaussian decomposition (purple lines) of the SH^ 
(1,2,5/2 <- 0,1,3/2) (top), 13CH+(1 ^ 0) (middle), and 
CH'*" (1 <— 0) (bottom) absorption spectra observed towards 
DR21(OH). The solid red lines and the dashed blue lines cor- 
respond to the confirmed (C) and uncertain (U) Gaussian com- 
ponents, respectively (see main text). The solid green lines are 
the empirical modelling of the CH"^ saturated line profiles. 



The results of the Gaussian decomposition procedure applied 
to the spectra are given in Tables lATI In turn, Figs. lA.ll - IA.SI dis- 
play the comparison between the result of the fit and the original 
data, along with the associated residuals. The uncertainties given 
in the Table I A. II are the formal 1 -cr errors derived from the di- 
agonal elements of the covariance matrix and do not take into 
account the systematic errors introduced by (1) the finite veloc- 
ity resolution, (2) the uncertainty on the continuum level, and 
(3) the error introduced by the decomposition procedure, which 
depends on its convergence criteria. 

According to i Godard et alJ (1201 Ol) (Appendix B), the finite 
velocity resolution of the spectra introduces an error on the 
linewidth determination smaller than 8 % for the CH^ and '^CH^ 
components, and smaller than 14 % for the SH^ components, if 
Al/ > 1.3 km s"' (the smallest observed linewidth). Furthermore, 
a relative uncertainty e on the continuum temperature leads to an 
error 6t ~ e-ln{l + ee^) on the opacity that ranges between 6% 
and 30 % when t varies between 0.03 and 2.7 (the lowest and 
highest observed central opacities) and using e - 10 % (value 



1.15 
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1.05 
1.00 
0.95 
0.90 

0.00 
-0.10 
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i ^ \ 




0.00 
-0.08 
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Figure A.2. Same as Fig. lA.ll for the observations towards 
G34.3H-0.1. 



corresponding to the error on the beam efficiency and the side 
band ratio, Roelfsema et al. 2012). 

Unlike the other numerical fitting method^^that are devised 
to find the global minimum of x^, the Levenberg Marquardt al- 
gorithm is a so-called single-shot method: used when the initial 
conditions are easily fixed, this algorithm takes advantage of the 
information on the derivatives with respect to each parameter 
to converge on the nearest local minimum. However, since the 
minimum only corresponds to a statistical estimate of the fit pa- 
rameters, the code can converge on different solutions close to 
each other depending on the initial condition, the x' valley to- 
pography, and the criteria to stop the optimization. To estimate 
the associated error, we analysed the dispersion of the best-fit 
solutions when varying the initial conditions: for each spectrum, 
we assumed an uncertainty of 2 km s"' on the initial guess of 
the positions and the linewidths of the Gaussians. The resulting 
4200 fits are found to bracket the optimum solutions given in 
Table lA.ll with a standard deviation of 10 % on the central opti- 
cal depths and the Gaussian linewidths. 

When combined, the robustness of the decomposition proce- 
dure, the finite resolution, and the uncertainty on the continuum 
temperature correspond to maximal errors on the calculation of 
Av and t of 20 % and 30 %, respectively. 



such as the grid-Powell method, the Montecarlo-Powell search or 
the simulated annealing method 
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Figure A.3. Same as Fig. lA.ll for the observations towards Figure A.4. Same as Fig. lA.ll for the observations towards 
W31C. W33A. 
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Figure A.5. Same as Fig. lA.ll for the observations towards 
W49N. 



Figure A.6. Same as Fig. lA.ll for the observations towards W5 1 . 
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Figure A.7. Same as Fig. lA.ll for the observations towards Figure A.8. Same as Fig. lA.ll for the observations towards 
SgrA*+50. SgrB2(N). 
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Table A.l. CH* (1 <- 0), "CH* (1 <- 0), and SH+ (1 <- 0) absorption line analysis results. 

CH+ (1 ^ 0) "CH+ (1 ^ 0) SH+ (1,2,5/2 ^ 0, 1,3/2) 

vo Al/ ± (t(Av) To + cr(To) i^o Au ± o-(Au) tq + cr(To) vo Av ± (t(Av) tq + o"(ro) W 

(kms"') (kms"') (kms"') (kms ') (kms"') (kms"') 



DR21(OH) 



-10.8 ± 0.3 


6.9 ± 0.5 


9.1 (-2) 


+ 


4 (-3) 


U 


-2.4 + 0.1 


3.7 + 0.1 


2.3 (-1) 


+ 6 (-3) C 


-2.6 + 0.1 


2.0 + 0.4 


1.8 (-1) 


+ 


3 (-2) C 


20.2 ± 1.8 


7.0 ± 2.4 


4.7 (-2) 


+ 


l(-2) 


U 


3.4 + 0.2 


5.1 + 0.4 


5.6 (-2) 


+ 4 (-3) C 


10.6 + 1.2 


5.2 + 2.7 


2.2 (-2) 


+ 


2 (-2) U 














10.1 ±0.4 


6.6 + 0.7 


8.2 (-2) 


+ 4 (-3) C 
















G34 


1.8 ±0.1 


5.0 ± 0.1 


4.6 (-1) 


± 


9 (-3) 


c 


25.0 ± 0.3 


4.8 ± 0.7 


5.1 (-2) 


+ 6 (-3) 


C 


27.5 


+ 0.3 


3.8 + 


0.6 


6.2 (-2) 


+ 


2 (-2) C 


9.2 ±0.1 


4.6 ± 0.2 


7.3 (-1) 


± 


l(-2) 


c 


31.3 ± 0.3 


3.5 ± 0.7 


4.2 (-2) 


± 7 (-3) 


C 


49.4 


+ 0.8 


7.3 + 


2.4 


3.3 (-2) 


+ 


2 (-2) U 


12.2 ± 0.3 


2.2 ± 0.9 


1.8 (-1) 


± 


2(-l) 


u 


38.8 ± 0.4 


1.7 ±0.8 


2.4 (-2) 


+ 1 (-2) 


U 


60.4 


+ 0.1 


2.2 + 


0.3 


1.1 (-1) 


+ 


3 (-2) C 


14.5 ± 0.5 


3.7 ± 1.6 


3.6 (-1) 


+ 


9 (-2) 


u 


42.1 + 0.2 


1.9 + 0.4 


5.2 (-2) 


+ 9 (-3) 


C 
















18.5 ± 0.6 


4.7 ± 0.9 


3.2 (-1) 


± 


4 (-2) 


c 


49.8 + 0.2 


6.2 + 0.4 


1.1 (-1) 


+ 6 (-3) 


C 
















24.8 + 0.2 


4.2 ± 0.5 


9.5 (-1) 


± 


5 (-2) 


c 


55.1 + 0.1 


3.5 + 0.3 


1.2 (-1) 


± 8 (-3) 


C 
















27.7 ± 0.2 


2.7 ± 0.4 


7.1 (-1) 


± 


2(-l) 


c 


59.8 ± 0.1 


2.8 + 0.2 


1.7 (-1) 


+ 9 (-3) 


c 
















30.2 ± 0.2 


4.0 ± 0.2 


1.3 (+0) 


± 


8 (-2) 


c 


62.9 + 0.1 


2.3 + 0.2 


9.9 (-2) 


+ 9 (-3) 


c 
















39.8 ±0.3 


4.8 ± 0.3 


6.0 (-1) 


± 


4 (-2) 


c 


























42.0 ± 0.0 


2.3 ± 0.2 


5.9 (-1) 


± 


8 (-2) 


c 


























61.3 ± 0.0 


1.9 ± 0.1 


7.1 (-1) 


± 


3 (-2) 


c 


























63.9 ± 0.1 


3.5 ± 0.1 


6.5 (-1) 


± 


1 (-2) 


c 


























W31C 


-9.9 ± 0.1 


5.0 ± 0.2 


1.4 (-1) 


± 


4 (-3) 


c 


-1.5 ±0.1 


3.5 + 0.2 


8.2 (-2) 


+ 5 (-3) 


c 


18.1 


+ 0.4 


4.2 + 


1.0 


2.7 (-2) 


+ 


7 (-3) C 


-1.8 ±0.8 


3.7 ± 0.6 


5.8 (-1) 


± 


3(-l) 


u 


8.3 ± 0.2 


5.9 + 0.4 


9.1 (-2) 


+ 4 (-3) 


c 


22.6 + 0.2 


2.8 + 0.4 


5.8 (-2) 


+ 


1 (-2) C 


0.0 ± 0.2 


2.8 ± 0.5 


7.2 (-1) 


± 


4(-l) 


u 


18.1 ± 0.2 


6.8 ± 0.5 


1.4 (-1) 


+ 5 (-3) 


c 


27.7 


+ 0.5 


4.6 + 


1.3 


4.5 (-2) 


+ 


8 (-3) C 


1.7 ± 0.1 


1.3 ± 0.2 


2.9 (-1) 


± 


8 (-2) 


u 


23.1 ± 0.4 


4.2 + 0.8 


3.4 (-2) 


+ 7 (-3) 


c 


32.9 


+ 0.4 


4.4 + 


1.3 


4.5 (-2) 


+ 


8 (-3) C 


55.2 ± 1.6 


4.9 ± 1.6 


7.6 (-2) 


± 


3 (-2) 


c 


28.0 ± 0.2 


5.7 ± 0.4 


1.5 (-1) 


+ 5 (-3) 


c 


38.3 


+ 0.4 


3.1 + 


1.0 


5.4 (-2) 


+ 


1 (-2) C 


57.8 ± 0.2 


3.1 ±0.4 


1.8 (-1) 


± 


6 (-2) 


c 


33.1 ± 0.1 


4.1 ± 0.2 


1.6 (-1) 


+ 7 (-3) 


c 


41.0 


+ 0.3 


2.2 + 


0.5 


5.1 (-2) 


+ 


3 (-2) U 














38.5 ±0.1 


3.9 ± 0.3 


4.2 (-1) 


+ 1 (-2) 


c 




























42.1 ± 0.2 


3.2 ± 0.6 


1.8 (-1) 


+ 2 (-2) 


c 




























46.0 ± 0.2 


3.7 ± 0.3 


1.7 (-1) 


± 9 (-3) 


c 
















W33A 


2.5 ± 0.1 


5.2 ±0.2 


2.8 (-1) 


± 


l(-2) 


c 


20.9 ± 0.8 


4.7 ± 1.3 


4.2 (-2) 


+ 1 (-2) 


c 


21.5 + 0.2 


2.4 + 0.5 


4.9 (-2) 


+ 


1 (-2) C 


10.1 ± 1.2 


8.6 ±4.9 


3.5 (-2) 


± 


5 (-3) 


u 


24.1 ± 0.1 


1.3 ±0.1 


1.4 (-1) 


± 1 (-2) 


c 


30.6 


+ 1.0 


2.7 + 


1.4 


8.6 (-2) 


+ 


7 (-2) U 


16.4 ± 0.5 


2.1 ±0.6 


2.4 (-1) 


± 


6 (-2) 


c 


28.9 ± 0.2 


5.0 ± 0.4 


1.9 (-1) 


± 7 (-3) 


c 


32.9 


+ 0.9 


2.4 + 


1.1 


8.7 (-2) 


+ 


9 (-2) U 


17.8 + 0.3 


1.5 ± 0.4 


1.4 (-1) 


+ 


1(-1) 


u 


33.2 ± 0.2 


2.7 ± 0.5 


1.4 (-1) 


± 1 (-2) 


c 


38.0 


+ 0.3 


2.1 + 


0.7 


3.6 (-2) 


+ 


2 (-2) U 














36.5 ± 0.2 


2.3 ± 0.6 


7.0 (-2) 


± 1 (-2) 


c 


41.6 


+ 0.2 


3.4 + 


0.3 


8.9 (-2) 


+ 


1 (-2) C 














40.9 + 0.1 


4.6 ± 0.3 


1.6 (-1) 


+ 6 (-3) 


c 
















W49N 


-3.4 + 0.1 


4.1 ±0.2 


1.9 (-1) 


± 


6 (-3) 


c 


9.9 + 0.1 


6.2 + 0.3 


1.3 (-1) 


+ 7 (-3) 


c 


10.0 + 1.0 


6.3 + 2.0 


1.8 (-2) 


+ 


2 (-2) U 


1.6 + 0.1 


2.8 ± 0.2 


2.0 (-1) 


± 


9 (-3) 


c 


13.0 + 0.1 


4.1 + 0.2 


2.3 (-1) 


+ 7 (-3) 


c 


12.8 + 0.1 


4.0 + 0.2 


2.2 (-1) 


+ 


1 (-2) C 


22.9 ± 0.8 


4.8 ± 0.7 


9.9 (-1) 


± 


8(-l) 


u 


16.5 + 0.1 


2.9 + 0.2 


8.8 (-2) 


+ 1 (-2) 


c 


16.1 + 0.2 


4.0 + 0.5 


7.2 (-2) 


+ 


2 (-2) C 


27.8 ± 0.5 


6.5 ± 4.1 


1.7 (+0) 


± 


3(-l) 


u 


23.3 ± 1.2 


9.9 + 2.3 


2.5 (-2) 


+ 3 (-3) 


u 


23.3 


+ 0.6 


5.8 + 


1.3 


2.6 (-2) 


+ 


7 (-3) U 


34.3 ± 0.3 


5.2 ± 1.6 


4.0 (+0) 


± 


4(-l) 


u 


34.5 ± 0.2 


5.1 + 0.4 


7.5 (-2) 


+ 7 (-3) 


c 


32.9 


+ 0.5 


5.4 + 


1.0 


2.7 (-2) 


+ 


7 (-3) C 


40.0 ± 0.4 


4.6 ± 1.0 


2.7 (+0) 


± 


4(-l) 


c 


39.5 ± 0.2 


4.1 + 0.4 


4.7 (-2) 


+ 4 (-3) 


c 


40.2 


+ 0.4 


5.5 + 


0.8 


3.1 (-2) 


+ 


7 (-3) C 


45.5 ± 0.3 


5.8 ± 0.9 


2.1 (+0) 


± 


9 (-2) 


c 


44.6 ± 0.9 


7.3 + 1.8 


5.2 (-2) 


+ 4 (-3) 


c 


51.5 


+ 0.3 


4.0 + 


0.6 


3.4 (-2) 


+ 


8 (-3) C 


51.1 ±0.2 


4.1 ± 0.2 


1.1 (+0) 


± 


K-l) 


c 


51.3 ± 0.4 


5.8 + 0.9 


2.6 (-2) 


+ 4 (-3) 


c 


61.6 


+ 0.2 


4.7 + 


0.3 


6.5 (-2) 


+ 


8 (-3) C 


74.7 ± 0.1 


4.1 ±0.2 


1.1 (-1) 


± 


6 (-3) 


c 


57.7 ± 0.3 


6.0 + 0.6 


7.2 (-2) 


+ 9 (-3) 


c 


68.6 


+ 0.8 


5.2+ 1.6 


1.5 (-2) 


+ 


7 (-3) U 














63.4 + 0.3 


5.7 + 0.5 


9.5 (-2) 


+ 5 (-3) 


c 




























68.4 + 0.2 


3.6 + 0.4 


5.5 (-2) 


+ 6 (-3) 


c 
















W51 


1.0 ±0.1 


2.3 ± 0.1 


1.8 (-1) 


± 


9 (-3) 


c 


41.1 ± 0.7 


7.2+ 1.6 


2.4 (-2) 


+ 5 (-3) U 


47.9 


+ 0.2 


4.5 + 


0.6 


5.1 (-2) 


+ 


1 (-2) C 


5.8 ±0.8 


3.7 ± 0.7 


5.0 (-1) 


± 


2(-l) 


c 


49.3 ± 0.1 


1.0 + 0.2 


5.9 (-2) 


+ 1 (-2) C 


55.5 


+ 0.2 


2.5 + 


0.4 


5.6 (-2) 


+ 


2 (-2) C 


8.3 ± 0.4 


3.2 ± 0.4 


6.7 (-1) 


± 


2(-l) 


c 


54.5 ± 0.1 


4.6 + 0.3 


1.3 (-1) 


+ 6 (-3) C 


64.8 + 0.1 


2.0 + 0.3 


6.4 (-2) 


+ 


2 (-2) C 


12.7 ±0.1 


3.7 ± 0.4 


2.4 (-1) 


± 


8 (-3) 


c 


64.9 ± 0.3 


3.7 + 0.4 


1.3 (-1) 


+ 2 (-2) C 
















17.8 ± 0.2 


3.3 ± 0.5 


1.1 (-1) 


± 


8 (-3) 


c 


69.0 + 0.6 


5.2 + 1.0 


8.9 (-2) 


+ 9 (-3) C 
















23.4 ± 0.1 


5.2 ±0.4 


3.6 (-1) 


± 


7 (-3) 


c 


























28.9 ± 1.2 


4.7 + 1.9 


3.3 (-2) 


+ 


l(-2) 


u 


























39.9 ± 0.1 


3.5 ± 0.2 


1.8 (-1) 


± 


7 (-3) 


c 


























48.1 ±0.7 


6.2 ± 0.7 


1.1 (+0) 


± 


2(-l) 


c 


























49.5 ±0.1 


2.4 ± 0.3 


8.8 (-1) 


± 


2(-l) 


c 


























54.4 ± 0.2 


5.8 ± 0.4 


2.0 (+0) 


± 


K-l) 


c 


























59.5 ± 0.1 


1.9 ± 0.2 


2.5 (-1) 


± 


3 (-2) 


c 


























70.3 ± 0.1 


2.1 ±0.4 


2.7 (-1) 


± 


1(-1) 


u 


























72.2 ± 0.5 


3.7 ± 0.6 


3.5 (-1) 


+ 


5 (-2) 


c 



























SgrA*+50 



-204.4 + 0.5 8.9+ 1.1 1.1 (-1) + 1 (-2) C -157.9 + 0.4 6.4 + 0.9 1.0 (-1) + 1 (-2) C -136.9 + 0.2 9.4 + 0.5 3.2 (-1) + 3 (-2) C 
" Validity of the Gaussian detection: C = confirmed if Av > 3cr(Au) and To > 2.5cr(To); U = uncertain otherwise 
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Table A.l. continued. 



CH+(1^0) "CH+(1^0) SH+ (1,2,5/2 ^ 0, 1,3/2) 



(kms-') 


Au ± cr(Au) 
(km s"') 


Tn + Cr(Tf\] 




V" 


(kms-') 


Av ± cr{Av) 
(kms-') 


To 


+ 


o-(to) 




(kms-') 


Av ± cr{Av) 
(km s-') 


Tq 


+ 


0"(To) 




1 01 n 


+ 04 


6.7 ± 1.0 


1.3 (-1) ± 2 ( 


-2~1 




1 47 4 + 

- IH- / It 


0.2 


2.1 ± 0.5 


1.3 (- 




It J 


f 9"! 
\ Z-) 




1 9fi n 

- 1 Z\J.\J 


+ n 1 

It U.J 


3.8 ± 0.6 


1 6 (- 






J 


z.) 




1 fiQ S 


+ fl 


20.1 ± 1.6 


7.3 (-1) ± 2 ( 




c 


1 Jo. J It 


0.9 


7.5 ± 1.5 


2.0 (- 


1) 


It z, 


( Z-) 


c 


1 14 


+ n 4 

It U.H- 


4.5 ± 0.9 


1 1 f- 


1) 




4 ( 


91 


c 


1 S'^ 
-1 JJ.J 


-t- n 
± U.J 


11.4 ± 0.8 


1 2 (+0) + 1 ( 


-21 


c 


-ijU.j + 


0.8 


7.3 ± 1.7 


2.0 (- 


1 


+ J 


f 9'» 


c 


1 nfi s 

-iUO.J 


+ n 4 

It U.T- 


3.6 ± 0.8 




1 




4 i 


91 


c 


147 Q 


+ n 1 


3.0 ± 0.2 


1 2 (+(f) + 8 ( 


-21 


c 


1 1 + 

1 1 J.U + 


2.7 22.7 ± 10.5 6.7 (- 


2 


+ 8 

+ o 




u 


07 n 

-y 1 .u 


+ n 1 

^ U.J 


2.4 ± 0.6 


1 2 


r 




ft t 
u 


91 


u 


117 4. 


+ 04 


5.8 ± 0.7 


1 8 C+0") + 1 ( 

J-.O \ i\JJ _!_ ± 1 


-11 




ion 7 + 


0.7 


6.9 ± 2.4 


7 4 T- 


2' 


It 


f 91 


u 


fi7 n 

o / .u 


+ n 

^ u.u 


4.8 ±1.4 


8 4 ^ 


2 




4 ( 


91 


u 


1 10 Q 


+ 

It O.J 


5.0 ± 1.2 


2 3 f+O') + 3 ( 


-11 


c 


fifi n + 

UU.U It 


1.1 


4.4 ± 2.5 


3.3 (- 


2' 


+ 9 
It z, 


f 9'! 
( z.) 


u 


'54 4 


+ n 

It U.J 


4.5 ± 1.1 


9 1 (- 


2 




4 ( 


91 
z.) 


u 


ins 1 


+ 08 
It w.o 


5.5 ± 1.8 


1.4 (+0) ± 4 ( 


-11 


c 


S^ 7 + 
jj. / It 


0.7 


8.8 ± 1.6 


1 %(- 

I .o \ 


2) 


It 1 


( z.) 


c 


^S 9 
jj.z 


+ n 1 

It U.J 


3.0 ± 0.7 


1 1 f- 


1) 




J 


91 


u 


08 


+ S 

It U.J 


3.3 ± 1.1 


7.0 (-1) ± 2 ( 


-11 


c 


^9 + 

JZ,.U It 


0.6 


9.0 ±1.4 


8 9 


2) 


It i 


I z.) 


c 


94 9 


+ n 1 

It U.J 


2.6 ± 0.8 


9 5 f- 


2) 




J 


91 
z.) 


u 


Ofi 9 


+ 9 1 

It .i. 1 


19.5 ± 3.1 


1.3 (+0) ± 1 ( 


-11 


c 


10 4 + 

1 7.^+ It 


0.5 


4.7 ± 1.0 


8.3 (- 


2) 


+ 9 
It z. 


I z.) 


c 


7 
I .y 


+ 1 n 

It i .u 


9.0 ± 2.2 


6.3 (- 


2) 




-t 1 
J 


91 


u 


78 7 


+ 04 

It U.H- 


5.2 ± 1.3 


1.8 (-1) ± 5 ( 


-21 


c 


4 7 + 
-H. / It 


0.3 


6.3 ± 0.7 


1.8 (- 


1) 


+ 9 
It z. 


I z.) 


c 


1 ft 4 


+ 9 1 
It z. i 


11.7 ± 4.2 


8.3 (- 


2) 




-t 1 
J 


91 


u 


68 7 
uo. / 


+ 04 

It U.H- 


14.0 ± 1.3 


4.6 (-1) ± 2 ( 


-21 


c 


4 + 
v.H It 


2.5 


10.4 ± 4.9 


9.1(- 


2) 


+ 9 
It z. 


( 9"! 
( z.) 


u 


91 


+ n s 

It U.J 


4.5 ± 1.3 


1 6 f- 


1) 




i 
y 


91 
z.) 


u 


S^ 4 


+ 01 

It U. 1 


6.8 ± 0.2 


2 3 C+O^ + 9 ( 


-21 


c 


14 9 + 
it.ii It 


0.4 


2.1 ± 1.1 


8.1(- 


2' 


+ 4 
It t 


( 91 


u 


11 8 

J 1 .o 


+ n s 

It U.J 


8.9 ± 1.7 


2 4 C- 


1 




4 ( 


91 


c 


-46.4 


± 0.2 


2.2 ± 0.6 


2.5 (-1) ± 7 ( 
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4( 


-2) 


c 
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c 
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c 


60.2 
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" Validity of the Gaussian detection: C = confirmed if Av > 3cr(Ai;) and tq > 2.5cr(To); U = uncertain otherwise 
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Appendix B: Determination of the excitation 
temperatures and of the subsequent column 
densities 

The column densities given in Table [3] are derived assuming 
a single excitation temperature Tex for all levels of a given 
molecule as 



Sttv; 



Q.. A.., L 



~hv„lkT,, 



17 



T dv 



(B.l) 



where vq, gu, gi and A„/ are the rest frequency, the upper and 
lower level degeneracies and the Einstein's coefficients of the 
observed transition, Q{T^y) is the partition funtion, and c is the 
speed of light. 

To derive the CH^ and SH^ 7 = 1 «— line excita- 
tion temperatures Tex, we adopt the recent computation of the 
CH^ de-excitation rate coefficients by collision with He and 
jHam mami et al. 2009; Lim et al. 1999), and assume a generic 
SH^ de-excitation rate coefficient by collision with H and H2 of 
10"'° cm"* s"'. Because the corresponding critical densities at 
100 K, ncrit,He(CH+) ~ 6 X 10^ cm-^ «erit,.-(CH+) ~ 4 x 10^ 
cm"^, and nci-it,H^(SH^) ~ 7 x 10^ cm"-', are considerably higher 
than those of the diffuse ISM, collisions are irrelevant in the 
excitation of CH^ and SH^, all the more so collisions of CH^ 
with H2 are its destruction pathway. Therefore, the level pop- 
ulations of these two ions only result from radiative excitation 
and possible additional excitation processes during their chemi- 
cal formation. Neglecting the latter, we derive the CH^ and SH^ 
7 = 1 «— line excitation temperatures Tex, from statistical 
equilibrium, using the submillimetre interstellar radiation field 
measured at / = 45° and b - Q° with the Cosmic Background 
Explorer (iReach et al.lll995i) . We obtain low excitation temper- 
atures rex(835GHz) = 4.4 K, and rex(526GHz) = 3.0 K that 
agree with the lack of CH+ (2^1) and SH+ (21 ^ 11), 
(22 <— 1 1), (21 <— 10), and (22 <— 12) detection in absorption at 
the velocities of the diffuse medium along the same lines of sight 
(PRISMAS and HEXOS observations, to be published). Note 
that an error of about 1 K corresponds to an error on the CH^ 
and SH^ column densities smaller than 1 %. While these values 
of Tex are valid for the velocity components associated with the 
diffuse interstellar gas, they clearly correspond to lower limits 
for those associated with the SFRs molecular environments. In 
the latter case, since the PRISMAS lines of sight target bright 
IR and submm emitting sources, a complete description of the 
IR radiative transfer accross the molecular gas surrounding the 
sources is needed (De Luca et al., in prep.; Vastel et al., in prep.). 

For excitation temperatures of '^ '^CH^ and SH^ of 4.4 K 
and 3.0 K, respectively, Eq. lB.ll finallv becomes 



A^(CH+) = 3.11 X 10'^ J" Tdv 

7^ 



-2 



and 



) = 3.05 X 10'^ Tdv cm 



,-2 



A^(SH+) = 1.42 X 10'' Tdv cm 



,-2 



(B.2) 



(B.3) 



(B.4) 



10'^ cm ^ and a 



using the CH+ 7=1^0, i^CH+ 7 = 1 ^ 0, and the SH 
N,J,F- 1, 2, 5/2 «— 0, 1, 3/2 transitions, respectively. 

The CH+ line opacity for N(CU^) 
linewidth of 4 km s ' is of the order of 15 and the expected line 
emission is about 4 mK (below the noise level, see Tables [T]and 
|2]l. The intensities of the '■'CH^ and SH^ lines that are optically 
thin are even lower 



Appendix C: Tracers of molecular hydrogen in the 
diffuse ISM 

The methylidyne radical CH has often been used as a probe 
of molecular hydrogen because of its linear relation with H2 
observed in the local d i ffuse medium (Federman 1982) and 
in dark clouds ( Mattilal Il986h . Compiling the UV and vis- 
ible data obtained o n 48 line s of sight (Savaae et al. 19771 
Rachford et al. 2002"; Crane et al. 1995; Crawford 1995 ; lAUen 
1994- Gredeletal. ,1993.: Federman et al. 1994; Dan ks et al. 
i 19841; [Jenniskens et al.lll992 ). Liszt & Lucas (2002) de rived a 
mean ratio N(CH)/N(li2) ~ 4.3 x 10"^ while Sheffer et al.1 
(12008!) found A^(CH)/A^(H2) ~ 3.5 x IQ-^, using their obser- 
vations of the absorption profiles of H2 and CH on a sample of 
90 diffuse clouds. Although the latter could be used in this work 
to estimate A^(H2), based on the Herschel observations of the six 
hyperfine components of the CH 7 = 1/2-3/2 absorption line 
(G erin et al. 2010^, it raises several problems that have to be 
taken into account. 



• As quoted bv lLiszt & LucasI (l2002h and lSheffer et al.1 ( l2008l) . 

there is a significant scatter (a dispersion of about a factor of 
3) on the A^(CH)-A^(H2) relation observed in the local diffuse 
medium, up to A^(H2) =10^' cm"^. 

• Despite the similarity of the CH absorption profiles with 
those of other tracers of molecular hydrogen (PRISMAS ob- 
servations, to be published), a A^(CH)-A^(H2) correlation has 
never been directly me asured in the inner G alactic ISM. 

• ICrane et al.l (1 19951) and Lambert et al.l (11990 ) found that the 
linewidths of CH observed in absorption are highly vari- 
able from one diffuse cloud component to another They re- 
ported the detection of CH line profiles either similar to or 
more Gaussian and less broad than those of CH^. This result 
was then confirmed by iPan et al.] (l2005h . who divided their 
observed CH absorption features into two categories: those 
with CN-like and those with CH^-like line profiles. This im- 
plies that the production of CH in the interstellar medium can 
sometimes be linked to the production of CH^and depends 
only slightly on the molecular hydrogen abundance (Godard 
et al., in prep.). 

For these reasons, we derive the H2 column densities from 
observations of the 7 = 1 «— ground-state rotational transition 
of hydrogen fluoride (HF) performed with the Herschel/HIFI 
instrument. If HF is saturated or not available, H2 is deduced 
from CH, assuming a mean HF/CH abundance ratio of 0.4 
(see columns 5 and 6 of Tabl e IHl. Th e a bsorption Une analy- 
sis performed by N eufeld et a l. (2010b) and lSonnentrucker et al.l 
(2010) showed that (1) the HF profiles are similar to those of 
H2O, and (2) HF accounts for more than 30 % of the total 
fluorine nuclei in the gas-phase along these sight-lines. Since 
both results agree with the chemical models that predict that 
the production of HF does not depend on the dynamics of the 
gas and that HF is the main F-bearing species, we deduce the 
A^(HF) - A^(H2) relation from an analytical study of the fluorine 
chemistry. For the range of parameters defining the diffuse and 
transluscent gas (Snow & McCall 2006), 10 < < 5000 cm"^ 
< Av < 5, andn(CO)/[C] < 0.9, the abundance of HF is driven 
by only four chemical reactions (.Neufeld & Wolfir6i2009.) : 



F H- H2 

HF H- y 

HF + C+ 

CF+ + Q- 



HF H- H 
F + U 
CF+ + H 
F + C 



(ki), 
(&f), 
(^2), and 
(^3), 



(C.l) 
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where 

ki = 1.00 X 10-1° [exp(-450K/r) + 0.078 exp(-80K/r) 

+0.0155 exp(-10K/r)] cm^s-', 
^HF = 1.17 X 10-1° exp(-2.2Av) s-i, (C.2) 
k2 = 7.20 X 10-°'* (r/300/:)-° i5 cm^s-', and 
h = 5.20 X 10-°** (7/300/:)-° '*° cm^s-i 

are their respective reaction rate coefficients from the cal- 
culations and laboratory measurements dZhu et al. 2002[ 
Brown & Balint-KurtfcOOOl: iNeufeld et al .''2005: N ovotnv etdJ 
2005 h. Therefore, assuming that there is no depletion of F on in- 
terstellar dust, and that the carbon is fully ionized and the sole 
source of electrons, the HF and H2 abundances are related by 

»(HF) ^ [¥] 

«(H2) ( [C] k2 + (nF/nn)/ki + 1/2/h4i + ArzAs)' 

where [F] and [C] are the fluorine and carbon elemental gas 
phase abundances relative to H nuclei (assumed to equal those 
observed in the solar neighbourhood, ~ 1.8 x 10~^ and ~ 
1.4x l O-'irespectivelv.lSavage & Sembachll996llSofia & Meved 
l200lt ISnow et aP l2007h . wh is the gas densitv and /hj is the 
molecular fraction defined as /h, - 2«(H2)/nH- For «h > 30 
cm-^ and Ay > 0.2, because F reacts rapidly with H2, the left- 
hand side of the denominator of Eq. IC.3l can be neglected, more 
than 90 % of the interstellar fluorine being in HF and less than 
10 % in CF^. We thus assume a n(HF)/n(H2) abundance ratio of 
3.6x lO-l 

Although HF is preferred to CH here as a tracer of molecular 
hydrogen, relation |C3] has to be used cautiously because: (1) it 
has never been proven observationally, (2) this abundance ratio 
linearly depends on the fluorine element al abundance, which can 
vary with the Galactocentric distanc e on (feCunhal 120041 
iRudolph et alJ2006l:lAbia et alj2010l) . and (3) it also depends on 
the depletion of fluorine onto grains, which can become impor- 
tant in the molecular gas. We note that the H2 molecular fraction 
computed as 2N(H2) / Nh on the velocity intervals given in Table 
ID ranges between 0.05 and 1 with an average value of ~ 0.4 
(Neufeld et al., in prep.). While this value suggests that the dif- 
fuse medium sampled by CH^, '^CH^, and SH^ is on average 
weakly molecular, it does not preclude the existence within this 
gas of regions where /h^ = 1 . A complete analysis of HF and CH 
as tracers of molecular hydrogen will be performed by Godard 
et al. (in prep). 
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